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Ab^mct 

This  report  examines  the  electrical  properties  ice  in  the  frequency  range  of  0- 
107  Hz,  attempting  to  be  suitable  both  as  a  simple  and  clear  textbook  for 
students  and  non-specialists  and  os  a  comprehensive  review  of  recent 
developments  and  discoveries  in  the  field.  Corresponding  to  this  double  goal, 
the  report  consists  of  two  parts.  The  first  one  is  written  in  textbook  style  and 
contains  most  general  theoretical  and  experimental  results  essential  for 
understanding  of  unique  electrical  properties  of  ice.  The  theoretical  interpr^ation 
of  ice  conductivity  and  ice  dielectric  permittivity  is  based  on  ice  being  a  protonic 
semiconductor.  Jaccard's  elegant  model  is  iteed  to  mathematically  describe  the 
electrical  properties,  and  is  expanded  on  cases  of  ice  samples  having  finite  size, 
boundaries  and  interfaces,  and  an  inhomogeneous  eiectric  field.  The  statistics 
of  charge  carriers  in  pure  and  doped  ice  is  discussed  in  detail,  as  are 
experimental  techniques  for  measurements  of  conductivity  and  dielectric 
permittivity.  The  first  part  contains  a  comprehensive  review  of  experimental 
results  on  ice  conductivity,  ice  dielectric  permittivity,  mobility  and  eiectric 
charges  of  protonic  charge  oorriers  ond  activation  energies  of  their  generation 
and  motion .  The  second  part  includes  more  ccxnplicated  physical  phenomeno— 
proton  injection,  dielectric  crossover,  icefield  effect  transistor,  thermostimuiated 
currents,  theory  of  ice/metal  interface,  field  extraction  of  charge  carriers  and  a 
recombination  injection. 


For  conversion  of  St  metric  units  to  U.S./British  customary  units  of  meosurematt 
consult  ASTM  Standard  E3dO-89a,  SlanckJni  F^ti(»forUseofth0lntBmatlonal 
System  of  Unlls,  published  by  the  Anmricon  Society  ter  Testing  and  Materials, 
1916  Race  St.,  Philadelphia,  Pa.  19103. 


This  report  is  printed  on  paper  that  contains  a  minimum  of  50%  recycled 
material. 
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FOREWORD 


At  the  present  time,  thousands  and  thousands  of  people  around  the  world  deal  with  ice, 
snow  and  permafrost.  They  are  scientists,  educators,  engineers,  navigators,  meteorologists 
and  others.  While  a  small  fraction  of  these  people  contribute  to  the  knowledge  base  in  ice 
physics,  all  of  them  use  knowledge  from  it  frequently.  Moreover,  successful  applied  research 
is  based  upon  fundamental  science — one  more  reason  for  ice  specialists  to  have  a  textbook 
on  ice  physics  on  their  desks. 

The  first  modem  ice  physics  text  was  Fletcher's  book  on  Tlie  Chemical  Physics  of  Ice  (1970) . 
Fletcher's  book  is  in  typical  textbook  format;  it  is  reasonably  brief  and  easy  to  understand. 
He  touched  on  a  few  of  the  most  important  topics,  but  not  all  of  them. 

The  most  comprehensive  book  on  ice  physics  to  date  was  published  by  Hobbs  in  1974. 
Hobbs  considered  almost  all  of  the  basic  aspects  of  ice  as  understood  at  that  time.  Moreover, 
he  described  and  compared  several  (sometimes  opposing)  viewpoints.  Tliis  fimdamental 
and  rather  large  (837  pages)  book  is  commoitly  known  as  the  "Ice  Bible"  by  specialists  in  the 
field.  In  1974  and  1975,  two  monographs  on  ice  were  produced  by  John  Glen.  These  were 
briefly  and  clearly  written  and  reviewed  almost  all  ice-related  subjects.  This  work  was  (and 
in  some  respects  still  is)  a  magnificent  introduction  to  ice. 

Finally,  in  1981  Maeno  wrote  a  simple,  popular  book  for  the  express  purpose  of  attracting 
people's  attention  to  the  subject. 

During  the  past  20  years,  a  sigiuficant  amount  of  new  experimental  and  theoretical  work 
has  appeared,  dramatically  changing  our  views  on  ice  physics.  As  a  result,  we  are  now  able 
to  formulate  physical  laws  using  more  simple  and  direct  methods.  We  have  found  some  of 
the  physical  models  used  in  the  past  to  be  completely  wrong.  The  physics  of  ice  is  a  much 
better  developed  subject  than  it  was  20  years  ago. 

For  the  above  reasons,  we  feel  the  time  is  ripe  for  a  contemporary  book  on  ice  physics, 
incorporating  the  known  and  proven  with  almost  20  years  worth  of  material  not  covered  by 
previous  works. 

I  have  tried  to  prepare  a  "readable"  book,  and  not  one  that  requires  the  reader  to  be  a 
uiuquely  educated  person.  It  is  my  intent  to  present  the  material  in  such  a  way  that  any  reader 
attracted  by  the  title  Ice  Physics  will  be  able  to  comprehend  it.  This  is  quite  difficult  for  a  book 
dedicated,  not  to  a  particular  field  of  knowledge,  but  to  a  specific  material.  Indeed,  for  ice  it 
means  we  have  to  consider  a  wide  variety  of  subjects,  including  quantum  chemistry,  solid 
state  physics,  the  theory  of  elasticity,  iortic  conductivity,  synchrotron  x-ray  topography, 
crystal  growth,  the  physics  of  surfaces  and  more. 

The  primary  goal  is  to  produce  as  simple  a  book  as  possible  without  sacrificing  scientific 
accuracy.  Experimental  facts,  physical  ideas  and  theories  will  be  strongly  orgartized  and 
bound  together  cohesively.  The  reader  will  be  introduced  to  a  wide  variety  of  material  on  a 
step  by  step  basis.  IThen  the  picture  will  be  whole. 

To  accelerate  materials  publication,  this  book  will  appear  first  in  the  form  of  a  series  of 
joint  CRREL-Dartmouth  reports,  later  to  be  published  in  CRREL's  Monograph  series,  on; 

1.  The  structure  of  ordinary  ice 

Part  I:  "Ideal"  structure  of  ice.  Ice  crystal  lattice 
Part  11;  Real  structure  of  ice.  Defects 

2.  Electriccil  properties  of  ice 

Part  I:  Conductivity  and  dielectric  permittivity  of  ice 
Part  II;  Advanced  topics  and  new  physical  phenomena 

3.  Optical  properties 

4.  Electro-optical  effects  in  ice 


iii 


5.  Thermal  properti^ 

6.  Mechaiucal  properties  of  ice.  Elasticity  and  anelastic  relaxaticm.  Plastic  properties. 
Fracture  of  ice 

7.  Electromechanicidetiectsinice 

8.  Surface  of  ice 

9.  Other  forms  of  ice  and  their  properties 

10.  Ice  in  space 

11.  Ice  research  laboratories 

The  reports  will  be  prepared  in  a  sequence  convenient  to  the  author.  The  present  report 
is  the  first  in  the  series. 
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Electrical  Properties  of  Ice 


VICTOR  F.  PETRENKO 

PART  1.  CONDUCmTTY  AND 
DIELECTRIC  PERMITTIVITY  OF  ICE 


INTRODUCTION 

A  wide  variety  of  the  physical  properties  of  materials  is  generally  considered  to  be  their  “elec¬ 
trical  properties."  Among  them  are  dc  and  ac  conductivities,  dielectric  permittivity,  piezoelec¬ 
tricity,  electrostriction,  photoelectric  phenomena,  etc.  This  report  is  devoted  only  to  conductiv¬ 
ity  and  electrical  polarization  of  ice  in  frequency  ranges  from  0  to  10^  Hz.  The  rest  of  the  phe¬ 
nomena  associated  with  electric  fields  and  electric  currents  will  be  considered  in  future  reports 
in  this  series. 

We  will  begin  with  some  elementary  definitions  and  then  compare  the  electrical  properties 
of  ice  with  those  of  other  materials. 

According  to  the  definition  of  conductivity  a 

7=^£  (1) 

where  J  is  electric  current  density  and  £  is  electric  field  strength. 

Equation  1  expresses  the  well  known  Ohm's  law  for  isotropic  materials.  In  the  case  of  an 
anisotropic  conductor,  we  have  to  use  the  more  complicated  equation 

/i=Oij£,  (2) 

where  Ojj  is  a  second-rank  tensor  of  conductivity  and /i  and  £j  are  components  of  vectors  /  and 
£  (Ney  1985).  Here  in  eq  2  and  also  below  we  use  Einstein's  sirmmation  convention:  when  a 
letter  suffix  occurs  twice  in  the  same  term,  suiiunation  with  respect  to  that  suffix  is  to  be  auto¬ 
matically  understood.  As  we  wall  see  later,  ice  Ij,  demonstrates  almost  perfect  isotropic  conduc¬ 
tivity  despite  its  hexagonal  symmetry.  That  is  why  we  will  use  eq  1  more  often  than  eq  2. 

When  one  kind  of  electric  charge  carrier  wdlh  electric  charge  q  and  a  concentration  of  n  is 
present 

}  =qn  <VA>  (3) 

where  <  v^>  is  a  mean  drift  velocity  of  the  charge  carriers  in  the  electric  field. 


If  the  electric  field  is  not  too  sttmg,  so  that  there  is  no  iivcrease  in  temperature  of  the  cckkIuc- 
tor  lattice  or  the  dtarge  carrrers,  ttie  mean  drift  velocity  <  v^>  turns  out  to  be  a  linear  function 
of  the  electric  field  strength  £ 

<yci>  =  ^  £  (4) 

where  the  constant  of  proportiorulity  p  is  called  mobility.  In  both  electron  and  ionic  conduc¬ 
tors,  p  turns  out  to  be  constant  in  a  wide  range  of  field  strengths,  but  for  different  reasons.  In 
electron  conductors  (metals  and  semiconductors),  the  conduction  electrons  are  moving  in  a 
chaotic  motion  with  very  high  velocities  (=  10^  m/ s  in  metals  and  =  10^  m/s  in  semiconductors 
at  room  temperature).  An  application  of  a  "weak"  electric  field  barely  changes  the  mean  abso¬ 
lute  value  of  the  velocity,  since  the  mean  drift  additions  to  the  velocity  are  very  small.  This  is 
why  the  average  time  of  relaxation  of  the  electron  impulse  Tp  does  not  change,  whether  it  is 
defined  by  either  collisions  with  the  defects  of  the  crystcil  lattice  or  by  scattering  on  the  phonons. 
The  mean  drift  velocity,  acquired  by  electrons  in  this  process  between  the  collisions  in  the  di¬ 
rection  of  the  electric  field,  is 

<Vd>=—  XpO.E  (5) 

ntg 

that  is 

e  Tn 

p  -  — t:  =  const  (6) 

ntg 

as  long  as  Xp  and  are  independent  of  £ .  Comparing  the  formulas  1, 3  and  4  gives  an  obvi¬ 
ous  result 


a  =  l<jlpn.  (7) 

If  the  charge  is  transferred  by  several  carriers,  moving  independently  of  each  other,  then  their 
contributions  have  to  be  added  together 


_  K  _  _ 

/  =I  /i=£ 


1  =1 


l«?il  Pi  «i 


£ 


(8) 


where  /  =  1, 2, ...  fC,  fC  being  the  number  of  carrier  types,  and  the  resultant  conductivity  a 
K 

o  =  Z  Oi  (9) 

1=1 

comes  out  to  be  equal  to  the  sum  of  the  partial  conductivities  of  all  the  carriers  of  the  charge  Oj. 

In  the  ionic  conductors,  mobile  ions  jump  between  the  closest  stable  positions  in  the  lattice. 
A  simple  case  of  one-dimensioruil  motion  of  an  ion  is  shown  in  Figure  1.  To  move  one  inter¬ 
atomic  distance  in  the  direction  of  the  field,  the  positive  ion  would  need  activation  energy 
-  qEall,  and  to  move  against  the  field,  +  qEall. 

The  frequency  of  ion  jumps  in  the  direction  of  the  field  is  expressed  by 


I  k^T  I 


(10) 
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Figure  1.  Potential  energy  of  an  ion  in  a 
crystal  lattice  Was  a function  of  coordinate 
X  with  fE  ^  0)  and  without  (E  =  0)  electric 
field  E;  a  is  an  inieratomic  distance. 


and  against  the  held 


^2  =  V  exp 


Urn  +  qaE/2\ 

k^T  I 


(11) 


where  v  is  the  frequency  of  ion  oscillations  in  a  potential  well  (in  equilibrium).  Since  every  jump 
moves  the  ion  over  the  distance  a,  it  is  apparent  that  the  mean  drift  velocity  of  the  ion  in  the 
direction  of  the  field  becomes 


<yci>  =  a 


(/l  “  fz)  =  2v  a  exp~^’"^kBTgJj^^\ 

\2kBTf 


since  usually 


(12) 


qEa  <<2kBT. 

We  can  substitute  the  argument  of  sh  instead  of  sft  in  eq  12 


<fd>  = 


(Y£fi)exp-^n 
\  kBTl 


JksT 


E. 


(13) 


(14) 


So,  the  ion  mobility  is 

^  _  /y.'??  I  (15a) 

UbTI 

and  is  independent  of  E.  usually  stands  for  the  movement  activation  energy. 

Also  notice  that 

v  exp"^*"  -  Vh .  (15b) 

This  represents  the  probability  of  jumps  or  the  mean  real  number  of  jumps  per  second.  For  cha¬ 
otic  thermal  motion  of  this  kind,  the  Einstein  relation,  connecting  mobility  p  and  the  diffusion 
coefficient  D  is  also  true 
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Figure  2.  Frequency  dependence  of 
conductivity  of  pure  monocrystal¬ 
line  ice;  a  is  measured  in  £2'^  m~^, 
f  in  Hz;  T  = -10  °C. 


Figure  3.  Direct  current  density  I 
through  the  same  ice  sample  as  for 
Figure  2  as  a  function  of  time  t;  E  = 
Ofor  t  <  0,  E  =  const  Jbrti0;l  =  Efcr, 
+  a^expl-tl  Tpl).  Thed^nition  of  Tp 
will  be  given  below. 


D=^H.  (16) 

In  ordinary  electronic  and  ionic  conductors,  a  is  independent  of  frequency  up  to  very  high  fre¬ 
quencies  (while  ft) «  i,  where  x  is  either  the  time  of  relaxation  or  the  period  of  ion  vibration). 

However,  ice  behaves  in  a  much  different  way.  A  characteristic  frequency  dependence  of 
electric  conductivity  of  pure  monocrystalline  ice  at  -10°C  is  shown  in  Figure  2.  The  most  no¬ 
ticeable  thing  is  that  there  are  two  major  values  of  electric  conductivity:  low-frequency  Oj  and 
high-frequency  which  differ  by  three  orders  of  magnitude.  Dependent  on  temperature,  the 
characteristic  frequency  of  transition  cod  =  which  is  called  the  Debye  frequency  and  changes 

from  10^  Hz  (0°C)  to  zero  at  the  low  temperature.  (At  77  K  this  frequency  is  already  immeasur¬ 
able  in  pure  ice.)  One  of  the  outcomes  of  the  frequency  dependence  in  Figure  2  is  the  unusual 
time  dependence  of  the  electric  current  through  ice,  if  a  steady  electric  field  is  applied  to  it.  This 
dependence,  created  under  the  same  conditions  as  in  Figure  2,  is  shown  in  Figure  3. 

It  is  interesting  that  the  temperature  dependence  of  electric  conductivities  of  pure  ice  (Cj  and 
a„„)  can  be  closely  approximated  by  exponential  dependencies  such  as 


ooc  exp 


Ea 

2k  zT 


(17) 


where  the  conductivity  activation  energies  are  equal  to  about  1  eV  for  and  1.16  eV  for 
(see  the  Review  of  Experimental  Results  on  Ice  Conductivity  and  Dielectric  Permittivity  section). 

Such  relationships  of  electrical  conductivities  are  typical  for  electron  semiconductors.  Even 
the  absolute  values  of  Eg  are  very  close  to  the  value  for  silicon  (1 .14  eV  at  T = 300  K  [Long  1968]). 
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Also,  the  absolute  value  of  of  pure  ice  at  T  =  O^C  is  typical  for  many  pure  electron  semicon¬ 
ductors.  One  of  our  first  goals  will  be  explaining  these  observations. 

The  second  mjqor  electric  characteristic  of  ice,  its  polarizability,  is  also  unusual.  The  state  ^f 
the  dielectric  medium  in  the  external  electric  field  £  is  described  by  a  vector  of  polarization  P 

P=lP,  (18) 

where  the  sumnwtion  is  done  over  all  elementary  electric  dipoles  p;  present  in  a  unit  volume. 
In  other  words,  P  is  the  dipole  momentum  per  unit  volume.  The  following  equation  is  also  true 

P  =  I  ^  i  n  (19) 

where  and  ri  are  respectively  the  charges  and  coordinates  of  elementary  charges  present  in  a 
unit  volume.  In  many  cases  eq  18  is  more  useful  than  eq  19. 

In  nonferroelectric  crystals 

P=£oX£  (20) 

where  x  is  the  electric  susceptibility  of  a  material  and  is  the  dielectric  permittivity  of  vacuum. 
X  is  the  major  value  determining  the  ability  of  a  dielectric  medium  to  polarize.  More  frequently 
found  in  practical  use  is  the  relative  dielectric  permittivity  e 

z  =  l +Xi  D  =eo£p  (21) 

where  D  is  a  vector  of  electric  displacement. 

Polarization  of  water  is  determined  by  two  factors.  The  first  is  the  large,  constant,  elemen¬ 
tary  dipole  of  a  water  molecule.  A  part  of  these  dipoles  is  oriented  along  the  external  field  £  , 
causing  the  creation  of  a  sizable  polarization.  The  larger  the  initial  elementary  dipole  pi  and  the 
higher  its  concentration  in  the  medium,  the  larger  the  value  of  e  that  can  be  achieved. 

Secondly,  the  value  of  p,  is  increased  somewhat  because  of  rearrangement  of  the  charges  in¬ 
side  the  H2O  molecule.  However,  for  water  and  ice  this  is  the  less  sigruficant  mechanism  for 
creating  e  of  the  order  of  3  only.  Figure  4  shows  the  experimental  temperature  dependence  of 


Figure  4.  Temperature  dependence  of  relative  static  dielectric  permit¬ 
tivity  cf  ice  (Johari  and  Whalley  1981)  and  water  (Malmberg  and 
Maryott  1956).  Extrapolation  equations  found  in  these  papers  were  used 
to  draw  the  plot.  In  ice,  below  =  200  K,  a  peculiarity  called  "crossover"  may 
occur. 
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the  static  dielectric  permittivity  e  for  ice  and  water  around  0®C.  The  most  noticeable  fact  is  that 
for  ice  £  =  10^,  even  larger  than  for  water  at  dus  temperature.  This  is,  of  course,  evidence  that 
even  in  the  "rigid"  crystal  structure  of  ice,  the  H2O  molecules  (and  their  pi)  have  the  ability  to 
turn,  like  the  molecules  in  a  liquid  or  a  gas.  The  second  important  fact  is  that,  despite  the  lower 
density  of  ice  (i.e.,  a  smaller  concentration  of  water  molecules  or  elementary  dipoles),  it  has  a 
larger  e  value  than  water.  Apparently,  this  cannot  be  understood  without  considering  intermo- 
lecular  interactions.  In  the  dispersion  curve  e  (Fig.  5)  for  pure  ice,  obtained  under  the  same  con¬ 
ditions  as  diagrammed  in  Figure  2,  we  again  notice  two  values;  e^,  static  or  low-frequency  per¬ 
mittivity  and  high-frequency  permittivity  £„..  nie  transition  from  Es  to  occurs  again  at  the 
same  frequeiKy /]>  e« = 32  is  amditioned  by  electronic  and  ionic  polarization  of  water  molecules, 
i.e.,  small  displacements  of  charges  within  such  a  "rigid"  body  as  water  molecules.  Large  static 
permittivity  £3  can  be  explained  only  in  terms  of  the  water  molecule's  considerable  reorienta¬ 
tion  in  space. 

This  short  introductory  review  of  the  electric  properties  of  ice  yields  a  general  idea  about  what 
kind  of  properties  the  physical  model  of  ice  conductivity  should  describe.  After  an  introduc¬ 
tion  of  this  model,  we  will  get  acquainted  in  detail  with  vast  experimental  results  and  conduct 
an  analysis  and  comparison. 


NATURE  OF  ELECTRIC  CHARGE  CARRIERS  IN  ICE 

To  get  acquainted  with  modem  theoretical  descriptions  of  the  electrical  properties  of  ice,  we 
will  discuss  a  question  concerning  the  nature  of  charge  carriers  in  ice. 

Any  physical  model  describing  the  electrical  properties  of  ice  must,  first  of  all,  be  based  upon 
knowledge  of  charge  carrier  type.  One  could  assiune  a  priori  that  charge  carriers  are  either  con¬ 
duction  electrons  (holes)  or  ions  (as  in  water).  As  suggested  by  Chesnakov  et  al.  (1987),  because 
of  the  fractional  electric  charge  of  ions  in  ice  (see  the  Other  Charge  Carrier  Transfers  and  Their  Ef¬ 
fective  Charges  section),  the  latter  can  accept  or  return  electrons,  changing  the  sign,  while  stay¬ 
ing  charged.  So,  these  peculiar  radicals  (ions  ±  electrons)  may  also  be  put  in  the  category  of 
possible  charge  carriers. 

The  dispute  about  the  nature  of  charge  carriers  in  ice  (electrons  vs  protons)  has  a  rather  long 
history.  Back  in  1957  Decroly  et  cil.  (1957)  carried  out  an  experiment  that  was  for  a  long  time 
considered  to  be  decisive  in  proving  the  proton  r\ature  of  ice  conductivity.  The  idea  of  the  ex¬ 
periment  was  simple  and  elegant.  If  m  ice  (as  in  v  ater)  the  app^earance  of  charge  carriers  is  caused 
by  thermoactivated  reaction  in  the  bulk 
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2H20^:?H30"^+oh- 


(22) 


then,  when  an  electric  current  is  passing,  one  would  expect  hydrogen  to  be  produced  at  die 
cathode 

4H30"^  +  4e-»2H2T  +  4H20  (23) 

and  oxygen  to  be  yielded  at  the  anode 

40H--4e^02t  +  2H20.  (24) 

According  to  eq  23  and  24,  in  the  case  of  ionic  conductivity,  the  passage  of  every  four  elec¬ 
trons  must  be  followed  by  the  yielding  of  two  hydrogen  molecules  and  one  oxygen  molecule. 
That  is,  Faraday's  law  of  electrolysis  must  be  fulfilled.  If  the  amount  of  generated  gas  comes  out 
to  be  less,  then  one  should  seek  other  charge  Ccirriers.  Such  an  experiment  using  direct  current 
requires  specuil  electrodes  to  be  designed.  The  prci>lem  is  that  a  stable  electron  exchange  be¬ 
tween  ice  and  metal  electrodes  is  usually  impossible.  Direct  current  through  an  ice  sample  with 
metallic  electrodes  ceases  after  a  short  time  (the  reasons  for  this  will  be  discussed  in  detail  in 
the  Charge  Exchange  at  Ice/Metal  Interfaces  section).  That  is  why  Decroly  et  al.  used  layers  of  fro¬ 
zen  HF  between  ice  and  metal.  Such  layers  can  easily  exchange  charges  with  both  ice  and  metal 
electrodes.  The  experiments  were  carried  out  at  -10®C  and  it  was  shown  that  exactly  the  same 
amotmt  of  hydrogen  was  generated,  as  follows  ft'om  eq  23. 

On  this  basis,  they  concluded  that  conduction  in  ice  was  protonic.  This  result  satisfied  the 
overwhelming  majority  of  scientists  and  was  not  re-examined  until  1990  (though  it  should  have 
been  done  in  the  1960s).  By  that  time  (1960s)  the  large  surface  conductivity  had  come  to  be  widely 
known  at  temperatures  above  -25‘’C.  This  conductivity  is  conditioned  by  the  presence  on  ice 
surfaces  of  a  so-called  quasi-liquid  layer.  This  layer  decreases  in  thickness  when  the  tempera¬ 
ture  decreases  (starting  from  ~1000  A  at  0°C),  and  its  contribution  to  conductivity  becomes 
negligible  at  T  <  -25  to  -30°C.  Besides,  if  we  take  into  account  the  extremely  great  mobility  of 
fluorine  ions  in  ice  reported  by  Haltenorth  and  Klinger  (1969),  we  have  to  agree  that  the  ice 
specimens  were  strongly  doped  and  most  of  the  impurities  were  being  contained  at  grain  bound¬ 
aries  and  surfaces,  thus  increasing  surface  conductivity  rather  than  bulk  conductivity. 

So,  the  conclusion  of  protonic  conductivity  can  be  made  only  with  respect  to  ice  surfaces  from 
the  experiments  of  Decroly  et  al.  Unfortunately,  even  this  inference  couldn't  be  drawn  from  their 
expreriments.  If  the  ice  specimen  ^'^ere  changed  to  a  piece  of  metal,  the  same  amount  of  hydro¬ 
gen  would  be  yielded  from  the ;  IF  layer,  which  is  a  protonic  conductor.  The  experiment  of 
Decroly  et  al.  proved  nothing,  although  indirect  evidence  of  the  protonic  conductivity  of  ice  ex¬ 
isted.  First,  the  dc  conductivity  of  heavy  ice  is  about  1(X)  times  less  than  that  of  ordinary  ice 
(Engelhardt  1973).  Second,  it  has  been  establidied  that  proton  injection  from  Pd  electrodes  into 
the  bulk  of  ice  increases  its  conductivity  (BuUemer  et  al.  1969,  Petrenko  et  al.  1983).  Certainly, 
these  results  do  not  mean  that  conductivity  is  entirely  determined  by  proton  motion. 

Eventually,  all  doubts  were  dispersed  after  Etecroly's  experiments  were  repeated  using  new, 
specially  designed  electrodes  that  provided  a  free  charge  exchange  between  ice  and  metal,  even 
at  very  low  temperatures  (Petrenko  and  Chesnakov  1990a).  We  w^ll  describe  in  detail  the  prin¬ 
ciple  of  action  of  such  electrodes  with  deposited  monomolecular  charged  layers  in  the  Charge 
Exchange  at  Ice/Metal  Interfaces  section. 

In  these  experiments,  hydrogen  and  ox)  ration  in  the  electrolysis  of  pure  ice  was 

measured.  The  dc  conductivity  of  pure  ice  was  entirely  determined  by  proton  motion  within 
the  temperatures  range  of  0  to  -42°C  (experimental  error  ±  3%). 
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JACCARiyS  MODEL  OF  ELECTRICAL 
PROPERTIES  OF  ICE  OF  INKNITE  SIZE 


At  present,  several  theories  describing  the  electrical  properties  of  ice  are  well  known.  We  will 
start  getting  acquainted  with  them  by  looking  at  Jaccard's  (1959, 1964)  model,  which  enjoys  a 
special  place  among  others,  for  several  reasons.  It  was  the  first  successful  model  of  ice  proper¬ 
ties  (both  dielectrical  and  conductivity)  and  is  extraordinarily  elegant  and  clear.  It  is  this  model 
that  is  disseminated  among  ice  physicists  at  the  present  time. 


Protonic  charge  carriers 

Jacccird  proceeded  from  the  assumption 
that  charge  carriers  in  ice  are  protons.  If  we 
take  a  look  inside  an  ice  lattice  (Fig.  6),  we 
will  discover  that,  for  protons  to  be  trans¬ 
ferred  at  sizable  distances  (such  motion  is 
necessary  for  dc  flow),  they  must  go  through 
a  succession  of  jumps,  both  along  hydrogen 
bonds  and  from  bond  to  bond.  It  is  obvious 
that  there  is  no  such  opportunity  for  protons 
in  a  defect-free  ice  structure.  Indeed,  any 
hopping  along  the  bond  from  one  molecule 
to  another  will  create  H30'*^  and  OH~  ions, 
and  any  jump  from  one  bond  to  another  will 
give  rise  to  a  pair  of  Bjerrum  defects:  L-  and 
D-defects.  Both  H30‘^  and  OH”  ions  and 
Bjerrum  defects  will  be  discussed  in  another 
report  (see  Foreword). 


Figure  6.  Chain  of  water 
molecules  in  the  ice  struc¬ 
ture.  Successive  hops  of  the 
proton  along  the  hydrogen 
bonds  and  between  them  are 
forming  charge  transfer  un¬ 
der  the  action  of  the  electric 
field. 


Since  the  creation  of  each  pair  of  such  defects  requires  energy  of  approximately  1  eV,  the 


process  of  proton  transfer  via  the  creation  of  defects  by  the  field  turns  out  to  be  possible  only  in 


extremely  strong  electric  fields. 


E  >  ieV  _  3  X  10^°  V/m  (25) 

e  a 

Actually,  in  a  real  situation,  the  application  of  an  electric  field  does  not  give  rise  to  new  protonic 
charge  carriers,  but  sets  in  motion  the  already  existing  ions  and  Bjerrum  defects.  The  structure 
of  hydrogen  bonds,  along  which  the  protons  move,  represents  itself  as  a  three-dimensional  grid. 
The  transfer  along  this  grid  of  any  of  the  above  four  charge  carriers  reorients  water  molecules 
in  such  a  fashion  that  the  second  defect  of  the  same  type  cannot  pass  along  the  same  path  in 
the  same  direction. 

To  understand  better  what  happens  to  water  molecules'  orientation  under  ionic  and  Bjerrum 
defect  transfer,  let  us  consider  the  track  of  an  arbitrary  charge  carrier,  a  positive  ion — H30'^  for 
instance.  When  we  are  speaking  about  ioruc  or  Bjerrum  defect  transfer,  it  is  important  to  keep 
in  mind  that  it  is  just  a  convention — as  a  matter  of  fact,  oxygen  atoms  do  not  hop  from  one 
molecular  site  to  the  next  one.  Only  protons  move;  plus  there  is  a  small  redistribution  of  elec¬ 
tronic  density  at  neighboring  molecules,  since  a  proton  is  a  powerful  center  of  attraction  for 
external  valence  electrons.  Thus,  when  speaking  about  H3O+  motion  we  have  in  mind  that  a 
characteristic  distribution  of  three  protons  near  the  oxygen  atom  shifts  in  space.  This  is  a  con¬ 
venient  convention. 

So,  let  us  make  an  ion  (H30'^)  track  into  a  straight  line  for  simplicity  and  follow  the  reorien¬ 
tation  of  water  molecules  as  the  ion  moves  from  left  to  right,  as  shown  in  Figure  7.  Comparison 
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Figure  7.  Schematic  representation  ofH^O^  ion  (a,  b)  and  D-defect  (c, 
d)  motion  along  a  chain  of  hydrogen  bonds  in  ice.  a)  Initial  arrangement 
of  H3O*  ion  and  water  molecules  in  the  chain  before  the  ion  moves  from  the  left 
to  the  right.  The  arrows  indicate  successive  proton  hops  (along  the  bonds)  by 
means  of  which  the  ion  moves,  b)  The  chain  after  the  ion  has  passed.  Obviously, 
no  one  H^O*  ion  can  go  along  the  same  chain  in  the  same  direction,  c)  D-defect 
motion  from  l^  to  right  along  the  same  chain.  The  arrows  indicate  successive 
hops  cf  protons  from  one  bond  to  the  next,  d)  The  chain  after  the  D-defect  has 
passed.  Comparison  with  Figure  7a  shows  that  the  motion  of  D-defect  restored 
positions  of  the  protons  along  the  chain.  It  is  ready  (or  "open")  now  for  the 
passage  of  another  H^O*  ion. 


of  Figures  7a  and  7b  leads  us  to  two  simple  conclusions.  First,  the  ion's  ability  to  jump  right  or 
left  is  determined  by  the  orientation  of  the  molecules  on  the  left  and  right  of  the  ion.  Let  us  des¬ 
ignate  the  concentration  of  water  molecules  that  enables  an  ion  to  jump  to  the  right  Wfight 
to  the  left  Secondly,  ionic  motion  to  the  right  switches  all  molecules  on  the  way  from  Mnght 

to  njgft.  Of  course,  the  reverse  holds  for  motion  to  the  left.  These  simple  observations,  as  it  turns 
out,  are  sufficient  to  derive  equations  that  describe  ion  fluxes  in  ice. 

Let  us  consider  the  direction  of  ion  motion  from  left  to  right  as  positive,  and  ice  to  be  com¬ 
posed  of  tracks  such  as  those  shown  in  Figure  7,  which  of  course  is  schematic.  It  is  obvious  that 
total  ionic  flux  to  the  right  h  will  be  proportional  to  1)  the  frequency  of  protons'  attempts  to  jump 
along  bond  v^,  2)  ionic  concentration  «i,  3)  the  difference  between  water  molecule  concentra¬ 
tion  (Krighf-Mieft)  arid  4)  the  length  of  jump,  i.e.,  the  distance  between  t'vo  water  molecules  r^ 
(0-0  distance). 

;i««l  (nright-  «left  )vh  ^00  •  (26) 

In  fact,  instead  of  Mnght  and  n\^f^  we  should  use  their  ratio  to  the  total  concentration  of  water 
molecules,  «h20  which  is  equal  to  the  probability  of  finding  a  molecule  in  either  the  bright  nr  in 
the  nieft  group. 


/l  *  "1  ("right  -  «left)vh  foo/"H20- 
Since,  at  first,  before  the  current  passage 


(27) 


"right  =  "left 


(28) 


and  the  difference  between  these  concentraticKis  is  just  equal  to  the  difference  in  the  number  of 
jumps  that  have  already  occurred  to  the  left  and  to  the  right  in  a  unit  volume 


"right  "left  —  ^ 


junrqx 

jumps 

to  the  left 

-I 

to  the  right 

since  beginning 

since  begirming 

of  motion  , 

of  motion  ^ 

(29) 


then  it  is  obvious  that 

t 

"right  —  "left  “  —  f  yi(f)df  (30) 

'■oo  ^ 

where  integration  is  held  over  t. 

In  the  initial  moment  of  time  /  =  0,  the  water  molecules  had  no  distinguishing  orientation. 
For  the  sake  of  notation  simplicity,  let  us  designate  the  integral  in  eq  30  as  Q 

f 

(31) 

0 

Then,  substituting  eq  31  and  30  into  27 

yi(f)«-"ivin(f). 


Using  Vh  from  eq  15a  and  15b,  derived  for  such  ion  jumps  {a  =  r^y  we  obtain 


h  (*  ) 


ei"H20  r2 


(32) 


where  denotes  the  effective  charge  carried  by  the  ions.  Since  n  H2O  we  get  the  ultimate 

result  ’’oo 


/i(f)  =  9£2 


(33) 


where  factor  <|)  =  §  fcgT  Toq.  A  dimensionless  term,  ^  cannot  be  found  from  such  a  simple  sche¬ 
matic  consideration. 

This  simple  analysis  reflects  most  of  the  significant  physical  aspects  of  ion  transfer  in  ice,  but 
does  not  take  into  account  the  complex,  three-dimensional  structural  topology  of  the  ice  crys¬ 
tal  lattice.  The  strict  analysis  is  extremely  difficult  and  time-consuming  and  has  not  yet  been 
carried  out.  Jaccard  was  the  first  to  give  a  theoretical  consideration  yielding  ^  in  1959  and  again 
in  1964.  The  results  obtained  by  Jaccard  (1964)  gave 

%  =  I6/V3  =  9.238 


regardless  of  defect  types. 


(34) 


In  the  previous  pi^>er,  Jaccard  (1959)  found  foat  for  ion  transfer 


^  =  32/(3V3)»6.16  (35) 

and  for  ^rrum  defects  transfer 

^  =  3V3  =  5.2.  (36) 

Currently,  we  prefer  to  use  the  results  obtained  by  Hubnxan  (1979a,  b)  frcHn  the  aiudysis  of  a 
large  niunber  of  experimental  data 

^  =  3.85±  5% .  (37) 

Careful  theoretical  treatment  conducted  by  Nagie  (1974)  showed  that  the  predicted  theoretical 
anisotropy  of  the  dielectric  properties  of  ice  should  be  of  the  order  of  IfT*,  if  any.  As  we  will 
see  later,  it  is  the  constant  ^  that  )rields  tihe  magrutude  of  foe  dielectric  permittivity,  so  we  will 
ccmsider  it  to  be  isotropic,  i.e.,  scalar,  and  extend  our  particular  result  (eq  33)  to  foree  dimen¬ 
sions. 


(38) 

where 

♦  =3.85  it  bT  Too  (39) 

and  Q, ,  which  we  will  call  a  configurational  vector  from  now  cm 

(40) 

0 

Thus,  we  have  discovered  that  ion  fluxes  (eq  38)  in  ice  depend  upon  a  "prehistory"  of  ion  mo¬ 
tion  (eq  40).  Let  us  see  what  else  it  leads  to. 

Electrical  properties — one  type  of  charge  carrier 

What  kind  of  conductivity  and  dielectric-permittivity  properties  should  ice  with  one  type  of 
charge  carrier  (H30^)  ions  exhibit? 

Let  us  start  with  ice  behavior  in  a  constant  electric  field.  Let  an  electric  field  £  =  £(,  be  applied 
to  ice  in  foe  initial  moment  of  time  t  =  0  and  foe  ice  be  impolarized  at  f  =  0.  Let  us  find  current 
density  /  as  a  function  of  time — foe  total  current  is  a  sum  of  currents  caused  by  foe  electric  field 
and  special  flux  from  eq  38 

' 

J  =e\h 

I  /l  =(ei  £-<!)«)  5L;  oi=einini.  (41) 

Cl 

t 

£i=J  ji  £  =  £o 

L  0 
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The  soluticm  of  this  simple  system  of  equaticxis  yields  an  exponentially  decaying  current 


}{t)=  lo  exp  {-t /to) 

(42) 

(43) 

/o  =  Oi  E©  • 

(44) 

Thus,  a  stable  current  does  not  exist  in  this  system.  Relaxation  time  Tq  is  called  Debye  relaxation 
time.  As  icms  shift  subject  to  E,  an  opposite  flux  arises,  caused  by  an  increasing  number  of  wa¬ 
ter  molecules  oriented  in  such  a  fashion  that  ion  motion  against  the  field  is  stimulated. 

In  a  state  of  equilibriiun  (/  =  «>),  these  two  fluxes  cancel  each  other  completely.  In  a  particu¬ 
lar  sense,  the  value  Q ,  a  configurational  vector,  is  similar  to  a  force  field  acting  upon  ions.  The 
dielectric  properties  of  ice  are  determined  by  rising  polarization  P 

P  =  lein;  ^  =  P  =  =7  (45) 

dt 

where  the  summation  is  over  all  mobile  charges.  Using  eq  42  and  initial  condition  P  (0)  =  0,  we 
see  that 

P(f)=|T{f)d»  =7oTD[l-e-‘/H-  (46) 

0 

The  stationary  magnitude  of  polarization  (f  =  «») 

P=7o'Cd=^£  =eoX£  •  (47) 

<t> 

In  an  external  electric  field  there  also  arises  polarization  that  is  not  associated  with  ion  displace¬ 
ments,  but  is  caused  by  a  redistribution  of  charges  inside  the  molecule 

Pe=eoXe£-  (48) 

The  usual  time  taken  to  produce  such  polarization  is  comparable  to  a  period  of  internal  molecular 
vibrations.  Hence  this  polarization  is  very  "fast."  The  total  polarization  Pt 

Pt  =Pe +P  =£o(Xe+X)£=£oXt£-  (49) 

The  entire  relative  dielectric  permittivity  of  ice  is  then 

es=l+Xt  =l+X+Xe=E«  +  X  (50) 

where  we  denote  the  "fast"  part  of  £5  as  e,„-  This  part  of  the  permittivity  has  no  connection  to 
the  motion  of  defects  and  it  determines  the  dielectric  properties  of  ice  at  "higji"  frequencies  when 
a)»l/To.  From  eq  50  and  47 

(es-eoo)  =  ei^/eo1>-  (51) 

One  can  easily  check  that  65  -  is  about  10^. 
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Other  charge  carrier  transfers  and  their  effective  charges 

Let  us  look  now  at  what  happens  when  other  protonic  charge  carriers  move.  First,  we  drould 
introduce  a  particular  formalism.  Let  1, 2, 3, 4  represent  H30'^  and  OH*  ions,  and  D-  and  L-de- 
fects  respectively.  Thus  C3,  for  instance,  denotes  the  electric  charge  of  a  D-defect  and  >12  stands 
for  OH~  ion  concentration.  Figure  7  shows  what  happens  in  a  water  molecule  chain  (discussed 
above)  if  an  H3O*  ion  or  I>defect  moves  along  it.  As  can  be  seen  from  this  sketch,  the  effect  of 
a  D-defect  is  quite  opposite  to  tihat  of  an  H3Q*’  ion.  Using  the  same  consideration  as  in  Figure  7, 
we  can  see  that  the  motion  of  an  L-defect  causes  the  same  reorientahon  of  hydrogen  bonds  in  a 
chain  as  in  the  case  of  H30‘*',  while  an  OH“  ion  causes  the  opposite  reorientation,  all  defects 
moving  in  the  same  direction.  Therefore,  in  the  case  when  all  four  charge  carriers  are  present, 
we  should  take  into  account  all  four  fluxes  affecting  the  formation  of  the  configuration  vector 
D .  Keep  in  mind  the  sign  change  mentioned  above. 

t 

■  (52) 

0 

If  we  introduce  =  (1,  -1,  -1, 1),  i  =  1, 2, 3, 4,  then 

dl  .  (53) 

Accordingly,  the  direction  of  defect  fluxes  caused  by  molecular  reorientation  is  different  for 

different  defects.  Having  conducted  the  same  consideration  as  in  the  case  of  H30'*’  ions,  we  will 
find 

/i  =U£-ni(t.fl).  (54) 

Equations  53  and  54  express  in  mathematical  terms  the  fact  that  a  flux  of  particular  defects  in 
ice  is  determined  not  only  by  the  magrutude  and  the  direction  of  the  applied  electrical  field,  but 
by  the  history  of  motion  of  all  charge  carriers  as  well! 

Let  us  touch  on  the  question  of  the  value  of  a  charge  carried  by  defects.  Since  H30'^  and  OH“ 
ions  and  L-  and  D-defects  are  created  and  annihilated  in  couples,  it  is  obvious  that 

ei=-e2ande3  =  -e4.  (55) 

If  we  place  an  extra  "foreign"  proton  into  an  otherwise  normal  ice  lattice,  we  will  inevitably 
obtain  two  carriers  at  once:  an  H3Q''  ion  and  a  D-defect,  as  shown  in  Figure  8.  Hence  the  sum 
of  those  charges  equals  a  proton  diarge  e 

e-i+e^  =  e.  (56) 

To  find  out  how  an  elementary  proton  charge  is  shared  between  an  ion  and  a  D-defect,  we  must 
consider  two  circumstances:  1)  a  fraction  of  the  intermolecular  distance  that  a  proton  travels,  in 
case  of  ions  and  Bjerrum  defect  transfer,  and  2)  displacement  of  electron  density,  when  a  charge 
carrier  is  moving  along  the  chain. 

The  most  strict  theoretical  analysis  of  this  problem  was  conducted  by  Scheiner  and  Nagle 
(1983),  who  found  that 

e3=(0.36±03)e .  (57a) 
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The  other  way  to  find  ei  and  is  experiment.  As  we  have  seen  (eq  51),  the  charge  of  the  pre¬ 
dominant  defect  is  incorporated  into  the  formula  for  e^.  Using  this  fact  and  analysis  of  experi¬ 
mental  data  of  measurements  of  Cj  for  specially  doped  ice  (in  sudi  fashion,  that  either  ions  or 
Bjerrum  defects  predominated),  Hubman  (1979a,  b)  has  found 

C3=(0.38±0.01)c  (57b) 

which  is  also  remarkably  close  to  the  calculatioi\s  of  Scheiner  and  Nagle.  Other  experimental 
results  for  ej  and  Cj,  obtained  during  recent  years,  are  found  later  in  Table  3.  These  results  are 
close  to  those  represented  in  eq  57a  and  57b. 

General  consideration  of  Jaccard's  model  with 
four  kinds  of  carriers 

Let  us  now  try  to  find  a  mathematical  description  of  a  and  e  as  functions  of  frequency  for  ice 
contciining  all  four  kinds  of  charge  carriers,  with  charges  e^,  mobilities  Hj  and  concentration  (i 
=  1, 2, 3, 4).  We  will  designate  through  Oj  their  partial  conductivities. 

Oi=lei|Hini.  (58) 

The  main  condition  restricting  the  universality  of  this  consideration  will  be  a  requirement  of 
uniformity  of  ice,  i.e.,  the  ice  has  no  boundaries  (is  infinite)  and  the  distribution  of  charge  carri¬ 
ers  is  uniform 


grad  Mi  =  0  .  (59) 

Equation  59  defines  the  absence  of  diffusicm  fluxes  of  particles  and  the  independence  of  all  fimc- 
tions  from  coordinates.  We  will  also  ccwisider  that  ice  is  electrically  neutral,  i.e.,  there  are  no  space 
charges. 

Suppose  an  oscillating  electric  field  is  applied  to  ice 

£  =  E  exp  {/(Of}  (60) 

where  co  =  the  cyclic  frequency  (o  =  2rc/ 
f  =  time 

/  =  fT 

=  complex  value,  as  usucd. 


The  total  density  of  electric  current  /  is  the  sum  of  the  drift  current  and  the  displacement 
current 


/  ~  /dr  /dis 


Jdis  =  ■ 

at 

The  drift  current  results  from  fluxes  of  all  particles 

~  4  _ 

/dr  =  51  A 

1=1 


(61) 

(62) 


(63) 


/i  =(ei£-Tii  i  =  l,2,3,4 


(64) 


(65) 


Tte  system  equations  (eq  61-65)  fully  defines  our  problem  for  eight  unknown  functions  f^is, 
/dr/  //  ft,  j  1, )  2, )  3,  i  4-  We  will  look  for  solutions  for  some  function  /  (to,  f)  as  usual  in  the  form 

f{(ti,t)  =  f  (to)  exp{ia)f}  166) 

where  /  (to)  is  an  amplitude  of  oscillations  variable  / .  Since  the  system  made  up  of  eq  61-65 
is  linecir,  after  substitution  of  eq  66  and  canceling  coefficients  exp  {tout},  we  obtain  a  rather  simple 
system  of  algebraic  equations.  The  procedure  of  finding  solutions  is  routine  and  we  present  it 
here  in  ultimate  form 


o(to) 


r  iti)TD(o«-Os) 

=  i-  =  Os - 5 


1  -l  tOTo 


(67) 


£{(0) 


=  £=o  + 


£s-£o. 

1  -i  oj-tD 


We  have  introduced  new  terms 


Ooo=S  Oi 
i=l 

Oi  +02  03  +04 


(68) 


(69) 


(70) 
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(71) 


g  —  € 

e'  (co)  =  Re  [e  (w)]  =  £„  +  — ^ — - — .  (74) 

1  +  ((OTd)^ 

We  are  most  interested,  of  course,  in  analysis  of  eq  73  and  74,  describing  the  real  parts  of  o  and 
e.  Now  we  will  proceed  to  this  analysis. 

Analysis  of  Jaccard's  results 

L'espite  the  complicated  form  of  cq  73  and  74  at  first  glance,  they  have  very  simple 
asymptotics,  namely 


Figure  9.  Ice  sample  connected  into 
measuring  circuit  (a)  and  its  electri¬ 
cal  equivalent  circuits  (b-d). 
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lim  o' (co)  =  <%;  limo' (w)  =  Ooo; 

(i>— >0^  CO — 


(75) 


lim  e' (<o)  =  Eg;  lim  e,  (to)  =  £„; 

(O-»0;  CO-»oo; 

Plots  of  the  functions  (/(to)  and  e'(oj)  are  depicted  in  Figures  2  and  5;  solid  lines  are  drawn  for 
the  values  of  parameters:  Xq  =  5.06  x  10"^  s,  <J„  =  1.9  x  10"^  m"^,  Oj  =  3  x  10^  m"^, 

e„  =  3.16,  Eg  =  100.8. 

The  transition  from  low-frequency  (s)  to  high-frequency  (»)  values  occurs  in  the  vicinity  of 
the  Debye  frequency /□  =  (0D/2it  =  1/2  tc  Tq  =  3.14  kl^  (Fig.  2  and  5). 

For  clarity  we  can  devise  an  equivalent  circuit  having  the  same  frequency  response  as  a  cube 
of  ice  with  sides  of  1  m  (Fig.  9). 

If  we  take  a  look  at  eq  69,  we  see  that  high-frequency  conductivity  a„  is  determined  by  the 
greatest  partial  conductivity  among  Oj.  Usually,  it  is  04  in  pure  ice  for  T  ^  -50°C,  i.e.,  conduc¬ 
tivity  of  L-defects.  The  fact  that  o„  incorporates  all  Oi  in  an  additive  fashion  implies  that  at  high 
frequencies  (<b»(i)d)  charge  carriers  move  independently  of  each  other.  One  can  imagme  this 
situation  as  a  parallel  connection,  as  shown  in  Figure  10a.  At  frequencies  much  smaller  than  cod 
and  for  direct  currents,  reverse  values  of  conductivities  rather  than  conductivities  are  summed 
up  (eq  70)  for  ions  Oj  +  02  and  Bjerrum  defects  O3 + O4.  This  is  common  to  series  circuits,  as  shown 
in  Figure  10b.  Such  conductivity  Oj  is  controlled  by  the  smaller  of  two  conductivities:  ionic  0;^^ 
=  01  +  02  or  Bjerrum  defects  Ob  =  03  +  O4.  In  pure  ice  at  relatively  high  temperatures 
(T  <  -  50°  C),  it  is  ionic  conductivity.  The  character  of  o^  is  attributable  to  the  fact  that  in  the  case 
of  dc  conductivity,  carriers'  fluxes  must  be  stationary.  Since  £  is  constant,  then  D  also  must  be 
constant,  as  can  be  seen  from  eq  64.  This  is  possible  only  under  conditions  of 

yi-/2~/3+/4  =  0  (77) 

from  which  (noting  that  =  ei  (/i  -  /y  and  /b  =  Cj  1/3  -j^])  we  obtain 

/ion  =  ii.  (78) 

ei  ^3 

i.e.,  ionic  and  Bjerrum  defect  currents  are  approximately  the  same.  Under  this  condition  (eq  78) 
all  hydrogen  bonds  restore  the  iiutial  position  from  time  to  time.  As  we  have  seen,  if  an  ionic 
current  "locks"  chains  of  hydrogen  bonds,  then  a  Bjerrum  defect  current  "unlocks"  them,  and 
vice  versa.  Using  an  equivalent  electrical  circuit  in  ice,  shown  in  Figure  9b,  we  can  easily  obtain 
a  time  dependence  of  current  through  the  ice  specimen  when  a  rectangular  voltage  pulse  is 
applied  (E  =  0  when  f  <  0;  £  =  constant  when  t  ^  0),  portrayed  in  Figure  3.  Of  course,  we  could 
obtain  the  same  dependence  from  the  system  of  made  up  eq  61-65  at  given  initial  conditions. 
In  the  Relaxation  Times  of  Electric  Polarization  and  Electric  Fields  in  Ice  section,  we  also  will  con¬ 
sider  effect  of  boundary  conditions  on  the  electric  relaxation  in  ice. 

Let  us  analyze  now  how  we  can  account  for  a  large,  static  dielectric  permittivity  of  ice  Eg.  In 
Jaccard's  model,  dipole  moments  of  water  molecules  do  not  appear!  This  is  an  advantage  since 
the  concept  of  molecular  dipole  moments  in  ice  has  proven  to  be  a  complicated  issue  (Nagle 
1979). 

The  polarization  P  is  defined  in  terms  of  the  spatial  displacements  of  charge  carriers,  as  can 
be  seen  from  eq  45.  Large  electrical  polarization  arises  during  a  time  interval  of  about  Xp  after 
an  electric  field  has  been  applied  and  is  determined  by  displacements  of  the  majority  charge 
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Figure  10.  Equivalent  electrical  cir¬ 
cuit  for  high-frequency  conductiv¬ 
ity  a„  (a)  and  for  static  or  loio- 
frequency  conductivity  (b). 


carriers  (usually  Bjerrum  defects)  under  the  action  of  the  electric  field.  This  displacement  grows 
until  a  force  resulting  from  Q  balances  the  electric  field  force.  The  balance  is  achieved  within  a 
characteristic  time  interval  Tq,  which,  as  can  be  seen  from  eq  71,  is  also  defined  in  terms  of  con¬ 
ductivity  of  majority  charge  carriers.  It  is  interesting  to  compare  Tq  with  the  well-known  Max¬ 
well  relaxation  which  is  the  case  when  the  only  restoring  force  is  the  electric  field  originat¬ 
ing  from  electrical  charges  building  up 

'Tm  =eoe<»/o=o-  (79) 

At  T  =  - 10°  C  Xvi  comes  out  to  be  about  30  times  smaller  than  Xp!  In  other  words,  £2  is  not  such 
an  efficient  restoring  force.  It  is  useful  to  estimate  an  average  displacement  of  majority  charge 
carriers  <x>,  building  up  the  polarization  of  ice 


P  =EoX£  =^4  <^>  «4 

(80) 

<^^_eo(es-l) 

(81) 

CTI4 

For  typical  values  of  E  ~  10^  V  m“^  and  n^  ~  10^^^  m  ^  (-10°C),  we  obtain 


<x>  ~  10  ^  m  =  10  A  (82) 

i.e.,  <x>  is  quite  microscopic.  Notice  also  that  <x>  «  =  10"^  m  is  the  average  distance 

between  L-defects.  Hence,  we  see  that  the  defects  move  independently. 

In  the  conclusion  of  this  section,  we  will  derive  formulas  for  so<alled  Cole-Coie  diagrams 
(Cole  and  Cole  1941),  which  are  frequently  used  for  experimental  data  analysis.  From  eq  68 

e((o)  =  e'((o)  +  ie''((o)  (83) 

where  e'((o)  is  given  by  eq  74  and 


£s-£oo 

e"(a))= - CO  -  Xn. 

1  -t-co^Xp 


(84) 
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It  is  easy  to  verily  that  e'  and  e"  are  connected  by  the  relationship 


e  - 


es-e. 


(85) 


which  can  be  represented  graphically  as  a  semidrde,  portrayed  in  Figure  11. 

If  experimental  data  closely  fit  the  Cole-Cole  diagram,  it  is  strong  evidence  in  favor  of  an 
existing  relaxation  process  with  only  one  relaxation  time. 


Figure  11.  Cole-Cole  diagram. 


ELECTRICAL  PROPERTIES  OF  ICE  OF  HNITE  SIZE 

In  previous  paragraphs,  we  have  conducted  a  theoretical  treatment  of  the  electrical  proper¬ 
ties  of  ice  that  has  no  boundaries  and  is  absolutely  uniform.  In  practice,  this  is  an  unattainable 
ideal,  and  attempts  to  compare  such  a  "uniform"  theory  with  experimental  results  can  raise 
doubts  about  the  credibility  of  this  theory,  although  this  theory  plays  no  role  in  that. 

For  electrical  measurements  to  be  carried  out,  ice  specimens  must  be  put  into  particular  elec¬ 
trical  circuits.  This  cannot  be  done  without  creating  some  kinds  of  interfaces,  such  as  ice/ dielec¬ 
tric,  ice/metal,  ice/semiconductor,  ice/protonic  conductor  or  even  ice/plasma.  The  presence 
of  these  interfaces  causes  considerable  change  in  the  simple  dispersions  a(co)  and  e((o),  depicted 
in  Figures  2  and  5.  However,  as  we  will  see  later,  extra  dispersions  appearing  in  these  experi¬ 
ments  allow  us  to  determine  the  concentrations  and  mobilities  of  numerous  charge  carriers  in 
ice.  In  practice,  it  turns  out  to  be  extremely  difficult  to  accomplish  an  efficient  electric  charge 
exchange  between  metallic  electrodes  (which  are  electronic  conductors)  and  ice,  which,  as  men¬ 
tioned  above,  is  a  protonic  conductor.  Ways  to  achieve  such  an  exchange  will  be  discussed  later. 

It  is  much  easier  to  suppress  such  an  exchange  completely,  by  placing  a  thin  dielectric  film 
between  ice  and  electrode.  The  thinner  the  film,  the  smaller  its  contribution  to  measurements. 
A  thin  oxidized  layer  of  AI2O},  i.e.,  sapphire,  at  the  surface  of  an  A1  electrode,  or  Si02  at  the 
surface  of  a  silicon  electrode,  can  successfully  play  the  role  of  such  a  thin  film.  These  electrodes 
are  called  perfectly  blocking,  and  the  electrical  properties  of  ice  with  such  electrodes  were  theo¬ 
retically,  exhaustively  corrsideredby  Petrerdco  and  Ryzhkin  (1984a). 

Experimentally,  e'(a))  is  determined  from  the  capacitance  C(a))  of  the  capacitor  in  which  the 
ice  specimen  to  be  investigated  acts  as  a  dielectric 

£'(o,)  =  ^i£M  (86) 
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where  we  consider  a  rectangular  ice  specimen  of  length  L,  directed  along  the  v  coordinate  and 
cross  section  s.  Two  perfectly  blocking  electrodes  have  the  same  area  s  and  are  located  at  v  =  0 
and  A  =  L.  The  problem  of  finding  e'(c»))  is  then  reduced  to  finding  C(a))  for  a  plane  capacitor. 
Let  us  suppose  that  an  ice  sp>ecimen  has  ideal  blocking  electrodes  that  arc  connected  to  an  ac 
circuit  so  that  a  surface  density  of  electric  charge  on  the  blocking  electrodes  caries  with  time 
according  to  the  harmonic  law  ocexp(-/ci)/l.  We  can  describe  this  charge  as  a  function  of  time  t 
and  coordinate  a  as 


Po(a,  t)  =  [5 (A)  -  6(.a-L)]  Osc  exp(-/cot/  (87) 

where  5(a)  is  a  delta  function  and  is  an  amplitude  of  the  surface  charge  on  the  electrodes. 
From  the  capacitance  determination 

C  (co)  =  Re  [sCscA<p  (w)j  (88) 

where  A(p(co)  is  the  amplitude  of  the  potenticil  difference. 

A<p(f )  =  ((o)exp  (-lojf) .  (8^)) 


Now,  to  find  A(p(co)  we  ha\’e  to  write  a  closed  system  of  equations  for  all  the  quantities  connected 
with  the  electnc  field  in  ice.  There  are  ten  such  quantities:  £(.v,f )  and  Q(x,t)  are,  respectively,  the 
projections  of  the  electric  field  strength  and  the  configuration  vector  onto  the  A-axis;  j;i(A,/)  are 
the  concentrations  of  the  defects  (/  =  1, 2, 3, 4),  and  j,(x,t)  are  the  projections  of  the  densities  of 
the  defect  fluxes  onto  the  A-axis.  This  problem  is  different  from  the  one  considered  above  in  two 
w'ays.  First,  by  the  presence  of  a  space  charge,  we  have  to  add  the  Poisson  equation 


a.r 


4 

=  X  +Po(x,t) 

I  =1 


(90) 


where  A;/j  =  (»;  -  (ijo),  ii,o  being  the  original  concentrations  of  the  defects  when  £  =  0. 

The  second  difference  from  the  homogeneous  problem  is  the  appearance  of  diffusion  terms 
in  the  equations  for  the  defect  fluxes 

/)=(£,£- 11.  (i  =  l,2,3,4).  (91) 

,_2  5a 

Adding  four  continuity  equati'  ns 

^  =  0  ,  i  =  1,  ..  4  (92) 

dt  5  a 

and  an  equation  for  the  configuration  vector,  we  obtain 

(93) 

i  =1 

Thus,  we  get  the  closed  system  of  ten  equations  (eq  90-93).  This  system  can  be  solved  in  a  lin¬ 
ear  approximation  when  I  CjAcp  I  «  k^T,  when  we  can  substitute  Oj  in  the  two  equations  (eq  91) 
by  making  Ojo  =  constant.  In  this  case,  by  employing  the  Fourier  transform  in  x,  we  shall  switch 
over  to  a  system  of  algebraic  equations  that  are  easily  solvable.  We  are  not  going  to  frighten  the 
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reader  with  this  procedure — we  will  only  show  the  results  obtained  in  that  way.  The  case  of 
I  CjAcp  I  »  kgT  will  be  analyzed  below  in  the  Field  Effect  Transistor  Made  of  Ice  section. 

Low-frequency  limit-screening  lengths 

With  a)-»0  an  electric  field  inside  the  ice  specimen  is  distributed  in  accordance  with  the  pres¬ 
ence  of  two  screening  lengths,  I/kj  and  I/K2 


E{x)  =  ^ 
2eo 


^  (exp"*^!^  +  exp“'^l(^  ^))  +  (exp”^^ 

^00  Cs 


+  exp-^<^(i--^)) 


In  the  case  when  Bjerrum  defects  <ire  majority  charge  carriers 


(94) 


Ki  = 


el  ("30-»-^40) 


1/2 


<2  = 


1/2 


(95) 


Figure  12  shows  the  electric  field  £  inside  the  ice  specimen.  At  the  screening  length  ic  j  \  which 
is  determined  by  the  concentration  of  majority  carriers  (ny^ + n^o)  and  coincides  with  the  well- 
known  expression  for  the  Debye  saeeiung  length  xj*,  the  field  drops  down  from  the  original 
value  asc/£o^  to  Osc/^o^-  As  we  have  seen,  eg = 10^,  i.  e.,  xj'  is  a  chm^cteristic  length  at  which 
Es  attains  a  steady-state  value.  At  large  x  the  field  drops  down  to  zero  because  of  screening  by 
minority  charge  carriers,  i.e.,  ions  with  a  larger  characteristic  screening  length  x^^  determined 
by  ion  concentrations  («io  ”20)- 

Figure  13  shows  an  "apparent"  static  ((i>->0)  dielectric  permittivity  of  ice  e'  as  a  function  of  a 
specimen  length  L.  First  of  all,  notice  that  the  only  value  of  the  permittivity  for  ice  of  thickness 
L  «  Xj^  is  E^  so  that  the  Debye  dispersion  (Fig.  5)  disappears.  Notice  also  rhat,  using  experi¬ 
mental  dependence  £'(0,  L),  we  can  determine  Xj  and  X2,  i.e.,  («io  +  n2o)  and  (”30  "4o)- 

Ignoring  the  dependence  of  e  on  sample  thickness  L  can  result  in  significant  errors  when  e  is 
measured  at  low  temperatures  for  which  x“^  and  L  are  comparable.  Since  within  surface  layers 


a.  Only  Bjerrum  defects  are  taken  into  account. 


b.  Bjerrum  defects  and  ions  are  considered. 


Figure  1 2.  Distribution  of  static  electric field  E 
inside  the  ice  sample;  Ej  =  E2 = 

(cfter  Petrenko  and  Ryzhkin  1984a). 
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FiguretS. Dependenceofa^ ,  ..rent static 
dielectric  permittivity  e'on  ice  thickness 
L.  L-dejects  and  HjO*  ions  attend  (after 
Petrenko  and  Ryzhkin  1984a). 


Figure  14.  Dependence  of  apparent  static  dielectric  permittivity  e' 
on  circular  frequency  O).  L-defects  and  Hfi*  ions  attend  (after  Petrenko 
and  Ryzhkin  1984a). 
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of  thickness  the  permittivity  is  es/e*  times  smaller  than  in  the  bulk,  it  is  necessary  to  take 
the  dependence  of  e  from  L  into  account  for  L  5  65/600  =  (05-1.0)  nun  at  T  =  150  K.  How¬ 

ever,  that  was  a  typical  sample  thickness  used  by  many  researchers  to  determine  the  critical 
temperature  T^.  in  the  Curie- Veiss  law 


1 

T-Tc' 


(96) 


Ignoring  6s(L)  dependence  could  also  be  responsible  for  wide  scattering  of  the  experimental 
value  of  Tf.  (see  the  Review  of  Experimental  Results  on  Ice  Conductivity  and  Dielectric  Permittivity 
section). 


Frequency  dependence  of  dielectric  permittivity 

As  was  shown  theoretically  by  Petrenko  and  Ryzhkin  (1984a),  on  ice/dielectric  boundaries 
new  dispersion  "steps"  appeared  on  the  6'((o)  plot.  While  the  main  Debye  dispersion  occurs  at 
cod  =  1/td,  where  Tq  is  a  characteristic  time  of  relaxation  of  majority  charge  carriers,  three  other 
characteristic  ftequendes  correspond  to  relaxation  times  Xi  for  three  other  types  of  charge  car¬ 
riers.  In  the  paper  referred  to,  an  explicit  solution  was  found  for  two  types  of  charge  carriers 
(for  instance,  L-defects  and  H30*'  ions).  The  frequerKy  dependence  for  this  case  is  shown  in 
Figure  14,  with  some  changes.  The  characteristic  parameters  indicated  on  Figure  14  are 

(Ol  =  D4K2I  (OD  =  ^{ai/ef  +  04 /el)  (97) 
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-1 


(98) 


£2 


L^s  kjL 


£l 


^oo  1-exp 
Es  k2L 


(99) 


Notice  that  the  four  equatiorAS  above  contain  only  four  ur\known  quantities,  lumely,  Dj,  D4,  n  j 
and  04,  since  Kj,  K2,  aj  and  04  are  the  functions  of  these  variables.  Therefore,  we  are  able  to  de¬ 
termine  all  of  them  from  an  experimental  curve  £(o)).  Owing  to  the  Einstein  relation  between 
Dj  and  Pj  (m  =  I  ej  I D;  / k^T),  we  find  m  as  well.  This  method  was  used  by  Zaretskii  et  ai.  (1987a) 
and  Zaretskii  (1991)  to  determine  the  mobilities  of  L-defiecIs  and  of  H30'*^  ions. 


CONCENTRATION  OF  CHARGE  CARRIERS 

So  far,  we  have  discussed  conductivity  provided  by  protonic  defect  motion  in  ice  (ions  and 
Bjerrum  defects),  unrelated  to  how  such  charge  carriers  can  appear  in  ice.  In  fact,  three  sources 
of  protonic  charge  carriers  in  ice  are  known.  Two  of  those  exist  in  a  state  of  thermal  equilibrium. 
They  are,  first,  thermcd  generation  of  carriers,  and  second,  atoms  and  molecules  of  impurities. 
The  third  source  is  various  types  of  external  excitation,  such  as  strong  electric  fields,  ultraviolet 
radiation  or  high-energy  particle  beams.  We  will  begin  our  discussion  of  charge  carrier  concen¬ 
tration  with  a  consideration  of  pure  ice  in  a  state  of  thermal  equilibrium. 

Intrinsic  charge  carriers 

Mobile  defects  of  a  protonic  subsystem  in  ice  can  be  created  by  thermal  fluctuations  in  a  crystal 
lattice.  As  in  the  nugority  of  solid  materials  (excluding  some  solid  gases  with  extremely  weak 
van  der  Waals  bonds),  file  typical  energy  required  for  an  intrinsic  point  defect  pair  to  be  cre¬ 
ated  is  about  1  eV.  As  we  will  see  later,  the  activation  energy  in  ice  is  about  1  eV  for  H30'*'  and 
OH“  ions  and  about  0.68  eV  is  required  for  the  creation  of  L-  and  EWefects.  These  values  largely 
exceed  the  average  energy  of  thermed  vibration,  k^T = 0.023  eV  (-10°C),  and  by  primary  inspec¬ 
tion,  the  appearance  of  a  large  number  of  these  defects  is  unjustified.  However,  it  is  common 
knowledge  that  the  state  of  thermal  equilibrium  at  craistant  temperature  and  voliame  is  attained 
at  a  minimum  of  free  energy  F  rather  than  at  a  minimum  of  internal  energy  U. 

F=U-TS.  (100) 

The  appearance  of  defects  increases  internal  energy  U,  of  course,  but  owing  to  an  increase  in 
entropy  S,  the  equilibrium  is  attained  at  some  firute  value  of  point  defects  concentration. 

We  can  illustrate  this  point  with  a  simple  model  example.  Let  us  consider  a  crystal  lattice 
consisting  of  N  identical  points  (atoms,  molecules  or  elementary  cells).  A  point  defect  with  ac¬ 
tivation  energy  Eg  could  appear  at  any  of  the  points  with  equal  probability.  Then,  the  part  of 
free  energy  ascribable  to  the  appearance  of  defects  will  equal 

Fd  =  Eon-TSc  (101) 

where  in  turn  is  the  cortfigurational  entropy  of  defects  in  the  lattice 


Sc=fcBlnW 


(102) 
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and  W  is  the  number  of  arrangements  of  n  identical  defects  over  tfiese  N  identical  sites 


W  _ . 

nl(N-ny. 

At  equilibrium 

^  =  0. 
dn 

Substituting  eq  101-103  into  104  and  using  the  Stirling  formula 

ln(N!)  =  N  InN-N  (N>>1) 
we  obtain 


(103) 

(104) 


n 

N-n 


Assuming  Eg  »  kgT,  N  »  n,  we  get 


n  =  N  exp 


(105) 


(106) 


To  be  more  precise,  the  appecirance  of  a  point  defect  not  only  changes  the  configurational  en¬ 
tropy,  but  also  affects  the  vibrational  entropy  because  of  a  slight  change  in  the  elastic  lattice  vi¬ 
bration  spectrum.  This  vibrational  entropy  change  -Sy  would  lead  to  an  additional  factor  of 
expfSy/lrB}  in  106/  but  as  Sy  is  so  small  the  factor  is  very  close  to  unity  and  so  will  be  omit¬ 
ted  below. 

Since  the  number  of  possible  protonic  configurations  in  ice  is  rather  large,  the  calculation  of 
Sf.  in  ice  is  more  complicated  than  in  the  above  example.  For  ice  with  ionic  defects,  such  a  cal¬ 
culation  was  conducted  by  Ryzhkin  (1985).  Considering  ice  with  ions,  he  first  noticed  that  each 
water  molecule  could  be  in  one  of  fourteen  states:  six  orientations  of  neutral  water  molecules 
that  correspond  to  the  possible  placement  of  two  close  protons  at  four  hydrogen  bonds  (4!/2!2! 
=  6);  four  orientations  of  OH”  ions  and  four  orientations  of  the  H30^  ion.  Let  the  number  of 
molecules  in  these  states  equal  N;  (/  =  1, 2,...14);  then  the  full  number  of  arrangements  of  N 
molecules  over  these  states  equals 

(107) 

Among  that  number  there  are  "right"  configurations  with  one  proton  on  each  bond  and  "wrong" 
ones  with  two  protons  or  without  any  protons.  The  probability  of  finding  a  right  bond  equals 
1/2.  Therefore,  for  the  number  of  right  configurations,  we  obtain 

14 

W  =(l/2)^N!/n  (108) 

i  =  1 

where  2N  is  the  number  of  bonds.  In  an  isotropic  case  we  have 

=N2...  =N6  =  (N-2ni)/6;  N7  =  N8...  =Ni4  =  «i/4  (109) 


N!/n  Ni!. 

I  =1 
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where  Kj  is  ttie  number  of  H30^  or  OH~  ions  (the  concentrations  are  equal  in  undoped  ice  when 
a  space  charge  is  absent ). 

The  entropy  of  the  ice  with  ions  equals 


^  =  In  W  =  -2«i  ln|- 

il-(A/  -2«i)ln 

2{N-2nJ 

kB  'A 

ll 

3N 

This  coincides  with  Pauling's  (1935)  result  W  =  when  «,  ->  0,  i.e.,  in  defect-free  ice.  Mini¬ 
mizing  F  with  the  entropy  given  by  eq  110,  we  obtain  for  the  equilibrium  concentration  of  ions 


N  -Irti 


(111) 


and  since  £3, »  Tk^T 

«i  a  2.  w  exp!- (112) 
3  1  Tk^Tl 

As  we  can  see,  the  specific  structure  of  ice  manifested  itself  only  in  appearance  in  eq  106  of  a 
factor  of  2/3  before  the  exponent.  The  factor  1/2  is  in  the  exponent  because  ions  exist  in  pairs. 

A  general  case  of  the  presence  of  ions  and  Bjerrum  defects  was  considered  also  by  Ryzhkin 
(1985,  unpublished).  Let  us  denote  D-defect  concentration  as  «[>  and  the  concentration  of  L- 
defects  as  «l-  W  D-  and  L-defects  are  formed  by  thermal  fluctuations,  i.e.,  in  pairs,  then  obvi¬ 
ously 

no  =  «L  •  (113) 

Using  a  similar  way  as  we  have  seen  above,  Ryzhkin  found  that  in  the  presence  of  both  types 
of  defects 


5/(^8^)  =  -2x  ln(x)  -(l-2z)  ln(2  (l-2z)/3) 

-4y  In(2y)- 2(1-2 y)ln(l-2y)  (114) 

where 

X  =  Tii  /N  and  y  =  no  /N .  (115) 

Minimizing  a  free  energy  of  ice  with  the  entropy  from  eq  1 14,  we  obtain  for  Bjerrum  defect  con¬ 
centration 


”D  =  «L=  ^  exp(-EaB/2  fcBl’)-  (116) 

Do  superionic  transitions  and  the  superionic  state  of  ice  exist? 

While  deriving  formulas,  accounting  for  H3CP'  and  OH“  ion  concentration  (eq  111  and  112) 
and  Bjerrum  defect  concentration  (eq  115  and  116),  we  considered  activation  energies  £aj  and 
EgB  to  be  constants.  In  1985,  Ryzhkin  noticed  that  this  is  not  always  the  case.  In  particular,  in 
the  case  of  large  defect  concentrations,  tire  activation  energy  is  no  longer  independent  of  defect 
concentration  that  results  in  a  phase  transition  under  certain  conditions.  The  main  idea  imder- 
lying  Ryzhkin's  theory  is  that  the  energy  of  formation  of  the  defect  pair,  ions  H30'^  and  OH", 
for  instance,  can  be  split  into  two  parts. 
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Figure  15.  Logarithm  cf  the  ratio  of  ion  concentra¬ 
tion  n  to  mol^lar  concentration  N  as  a  function  cf 
reciprocal  temperature  T.  The  steep  rises  correspond  to 
transition  into  the  superwnic  state  of  ice.  The  upper  curve 
is  for  Ej/E2  =  0.26,  the  lower  curve  for  E2/E2  =  0.58 
(after  Ryzhtdn  1985). 


E,/E2 


Figure  16.  Temperature  of  transition  cfice  into 
superionic  state  T^asa function  of  E2/E2  (af¬ 
ter  RyzMdn  1985). 


£a=£i+£2(n)  (117) 

where  £1  is  the  energy  of  fonning  an  ion  pair  at  distance  r^g,  i.e.,  at  two  neighboring  water 
molecules.  In  other  words,  £1  is  the  energy  required  for  the  proton  to  be  moved  along  the  hy¬ 
drogen  bond  from  one  neutral  water  molecule  to  another  neutral  one.  E^n)  is  the  energy  equal 
to  the  work  of  separation  of  formed  ion  pairs.  This  part  of  the  activation  energy  can  be  damp¬ 
ened  drastically  because  of  the  screerting  effect.  (Icms  are  screened  electrostatically  by  other 
mobile  charge  carriers.)  In  foct,  this  must  result  in  a  kind  of  "chain  reaction":  increasing  ion 
concentration  (by  means  of  gradually  increasing  temperature)  will  lower  Eiin)  and  hence  the 
entire  £3,  which  in  turn  causes  an  increase  in  n  and  a  decrease  in  £,.  This  relationship  of  E^  and 
n  becomes  essential  only  when  the  concentration  reaches  some  critical  value,  which  is  followed 
by  an  £3  jump  and  an  irrarease  of  n  by  orders  of  magnitude  (see  Fig.  15).  Unfortunately,  it  is 
extremely  difficult  to  calculate  the  value  of  £2  precisely,  because  of  the  uncertain  value  of  di¬ 
electric  permittivity  e„,  at  distances  about  interatomic  distance.  That  is  why  Ryzhkin  assumed 
that  1 5  S  3.2,  i.e.,  it  lies  between  values  of  e  for  vacuum  and  average  macroscopic  value.  The 
result  of  this  xmcertainty  of  e  is  that  the  theory  can  predict  only  qualitative  dependence  of  criti¬ 
cal  temperature  as  a  function  of  ratio  E1/E2,  as  diown  in  Figure  16.  Since  the  transition  indi¬ 
cated  above  does  not  occur  at  7=  273  K,  this  implies  that  Ei/£2  ^  026  and  e  >  253. 

The  sharp  change  in  concentration,  indicated  in  Figure  15,  can  be  interpreted  as  a  transition 
into  the  superionic  state,  since  ion  concentration  reaches  =  0.1  - 1,  i.e.,  it  becomes  comparable 
with  molecular  concentration.  Of  course,  in  this  case  the  ice  structure  itself  can  experience  quali¬ 
tative  changes,  resulting  in  a  disordering  of  the  oxygen  system  (quasi-liquid  state). 

Naturally,  all  considerations  presented  by  Ryzhkin  for  ions  also  hold  for  Bjerrum  defects. 
Under  usual  conditions  the  superionic  transformation  predicted  by  Ryzhkin  does  not  occur  up 
to  ice  melting  temperature.  However,  this  effect  could  be  significant  on  the  ice  surfoce  or  in  the 
vicinity  of  the  dislocation  core,  where  charge  carrier  concentration  can  exceed  the  bulk  ccmcen- 
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tratim  by  several  orders  of  nu^nitude,  becaiise  of  the  presm£  of  str(x\g  electric  fields  ai\d  elastic 
strain  fields  (Ossipyan  and  Petrenko  1988,  Ryzhkin  1992). 

Protonic  charge  carriers  introduced  by  doping 

Although  liquid  water  is  a  very  good  solvent  for  many  substances,  its  solid  phase — ice — 
dissolves  to  almost  nodung.  Most  of  the  impurities  dissolved  in  water  either  get  repelled  from 
ice  back  into  the  w  ater  or  precipitate  as  a  second  phase  inside  the  ice  bulk,  in  the  form  of  segre- 
gaticms  and  clusters  during  the  process  of  ice  growth.  There  are  some  exceptions,  such  as  a  few 
adds  (HF,  HCl),  ammonia  (NH3),  some  alkalies  (KOH,  NaOH)  and  their  derivatives  like  NH4F 
or  KCl.  These  substances,  first,  can  be  included  properly  into  the  ice  crystal  lattice  and,  second, 
they  drastically  change  the  protonic  carrier  concentration  in  the  ice. 

To  understand  how  these  impurities  affect  tfre  electrical  properties  of  ice,  let  us  consider 
probable  patterns  of  their  incorporation  into  the  ice  lattice  (see  Fig.  17).  Note  that  although  there 
now  exist  a  number  of  experimental  techruques  for  investigating  real  atomic  structure,  which 
allow  us  to  obtain  atomic  resolution  (transmission  electron  nucroscopy  of  high  resolution,  scan¬ 
ning  tunneling  microscopy,  atomic  force  microscopy),  the  real  structure  of  such  defects  in  ice 
has  not  been  studied  yet.  Patterns  such  as  those  represented  in  Figure  17  turn  out  to  be  just 
hypotheses  (assumptions)  based  upon  iorric  radii  analysis,  as  well  as  arralysis  of  energies  and 
lengths  of  atomic  bonds,  etc.  The  ionic  radius  of  F",  for  example,  is  rather  dose  to  that  of  Cr^, 
and  it  is  natural  to  assume  that  fluorine  substitutes  for  oxygen  successfully  in  an  ice  lattice.  As 
can  be  seen  from  Figure  17a,  since  the  HF  molecule  can  provide  only  one  proton  for  the  four 
nearest  hydrogen  bonds,  one  bond  appears  to  be  lacking  any  protons  at  all,  i.e.,  substitution  of 


a.  HF.  Top:  Ice  lattice  with 
incorporated  HF  molecule.  Bot¬ 
tom:  The  same  after  H3O*  ion 
and  L-defect  were  released  by 
means  of  sequential  hops  of  pro¬ 
tons  (along  the  hydrogen  bonds 
to  release  H3O*  ion  and  between 
the  bonds  to  release  L-defect). 


Figure  17.  Incorporation  of  impurities  into  the  ice  lattice. 
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b.  NHy  Top:  NH^  molecule 
incorporated  into  restructure. 
Bottom:  The  same  after  release  of 
OH~  ion  and  D-defect  from  the 
NHj  molecule. 


c.  KOH.  Top:  Diagram  to  indi¬ 
cate  how  KOH  dissolved  in  ice 
is  thought  to  introduce d^ects. 
The  K*  ion  is  interstitial,  and  the 
OH~  ion  substitutes  for  a  water 
molecule  leaving  one  bond  with¬ 
out  any  proton — an  L-defect. 
Bottom:  The  same  after  L-defect 
and  OH~  ion  have  gone  off  the 
KOH  molecule. 


Figure  17  (cont'd).  Incorporation  of  impurities  into  the  ice  lattice. 
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the  water  molecule  by  the  HF  molecule  results  in  the  appearance  of  an  L-defect  assiviated  with 
the  impurity  molecule.  The  binding  energy  of  such  an  L-defect  with  the  HF  molecule  v)b  viously 
is  not  large  (<  eV,  Jaccard  1959,  Camplin  et  al.  1978).  Since  the  number  of  possible  .urange- 
ments  of  the  l.-defect  in  an  ice  lattice  is  enormous,  practically  all  L-defects  leave  then  parent 
molecules. 

Moreover,  it  turns  out  that  a  single  proton  in  the  HF  molecule  is  not  boimd  firmly  and  can 
"jump  over"  along  a  hydrogen  bond  to  another  water  molecule,  because  of  the  thermoactivation 
process,  turning  the  water  molecule  into  an  H3O*  ion.  It  Ls  important  to  emphasize  the  chemis¬ 
try  of  HF.  To  a  physicist  the  natural  thing  would  be  for  an  extra  proton  to  attach  to  HF  forming 
an  OH“  ion,  but  chemically  HF  is  a  weak  acid  in  water  and  prefers  to  lose  a  proton. 

Thus,  HF  doping  of  ice  increases  the  L-defect  and  H30^  ion  concentrations  and  lowers  D- 
defect  and  OH~  ion  concentrations.  The  latter  occurs  because  the  equilibrium  concentrations  of 
D-defects  and  OH~  ions  in  pure  ice  are  determined  by  a  balance  between  the  rate  of  their 
thermoactivated  production  and  the  rate  of  recombination:  D-  with  L-defi  cts  and  H30^  with 
OH~  ions.  HF  doping  increases  L-defect  and  H30^  ion  concentrations,  so  that  the  rate  of 
recombination  increases,  which  results  in  decreasing  D-defect  and  OH“  ion  concentrations. 
Algebraically,  this  can  be  described  in  terms  of  the  products  of  concentrations 
«i  •  K2  and  M3  M4  being  constant  in  a  state  of  thermal  equilibrium,  irrespectic  e  of  doping 

Ml  -  mt  =  (Mi)2  =  i  N  ^  exp  -^^1  Mi<</V  (118) 

9  ksTi 

^3  •  ^^4  =  (”d)^  =  N  ^  exp  |-  ,  mo<<N  .  (119) 

I  k^Tl 

Another  molecule  that  can  successfully  be  incorporated  into  the  ice  lattice,  substituting  for 
water  molecules,  is  ammonia,  NH3  (Fig.  17b).  Since  NH3  donates  three  protons  (instead  of  two 
as  in  the  case  of  H2O)  for  four  hydrogen  bonds,  we  observe  two  protons  at  one  bond  at  the  same 
time,  i.e.,  the  D-defect.  As  a  consequence  of  the  thermoactivahon  pro-'ess,  this  defect  can  be  re¬ 
leased  and  take  part  in  the  conductivity  process.  Besides,  the  NH3  molecule  can  accept  a  fourth 
proton,  forming  an  OH"  ion  at  a  neighboring  water  molecule.  As  will  be  shown  later,  acti\  a- 
tion  energies  for  the  release  of  D-defects  and  OH"  ions  from  the  NI I3  m  !■  "ule  in  ice  are  quite 
sizable.  Therefore,  although  NH3  doping  of  ice  increases  D-defect  and  OH"  ion  concentrations 
(decreasing  L-defect  and  H30^  ion  concentrations  in  accordance  with  eq  1 1 8-1 19),  its  influence 
in  terms  of  conductivity  is  not  as  efficient  as  the  HF  doping  of  ice.  It  is  also  important  that  NH3 
doping  increase  the  concentration  of  less  mobile  carriers  (D-defects  and  OH”  ions),  as  will  be 
shown  in  the  Review  of  Experimental  Results  on  Ice  Conductivity  and  Dielectric  Permittiviti/  section. 

Lastly,  the  third,  very  important  kind  of  impurity  that  considerably  affects  the  conductivity 
of  ice  is  alkali  hydroxides,  among  which  the  most  significant  is  KOH.  The  most  plausible  mecha¬ 
nism  of  incorporation  of  KOH  into  the  ice  structure  is  shown  in  Figure  17c.  The  potassium  atom 
is  accommodated  in  interstices,  while  the  hydroxyl  group  OH  is  built  into  the  lattice.  As  can  be 
seen  from  Figure  17c,  this  results  in  the  appearance  of  L-defects  and  OH"  ions  linked  to  KOH 
impurity.  Because  of  thermoactivation  processes,  they  can  be  released  into  the  ice  bulk,  increas¬ 
ing  L-defect  and  OH"  ion  concentrations  and  suppressing  D-defects  and  H30'^  ions.  Because 
of  the  extremely  small  binding  energies  of  L-defects  and  OH"  ions  with  KOH  molecules  in  ice, 
these  impurities  become  efficient  donors  of  L-defects  and  OH"  ions  even  at  very  low  tempera¬ 
tures  (see  the  Review  of  Experimental  Results  on  Ice  Conductivihj  and  Dielectric  Permittivity  se^Hon). 

Let  us  briefly  discuss  the  question  concerning  the  statistics  of  charge  carriers  in  doj  cd  ice. 
This  problem  has  been  treated  many  times  in  scientific  literature  (see  Jaccard  1959,  Kroger  1974, 
Camplin  et  al.  1978).  Unfortunately,  in  all  instances  known  so  far,  this  consideration  was  based 
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upon  the  law  of  mass  action,  the  specifics  of  ice  being  disregarded,  which  can  yield  only  ap¬ 
proximate  results.  The  dependence  of  defect  conoentratioi\s  on  doping  impurity  ccsvcentraticms 
in  ice  was  considered  in  the  greatest  detail  by  Kroger  (1974).  We  will  give  below  an  example  of 
such  an  aixalysis  conducted  for  ice  doped  witi\  HF.  Substitution  of  water  molecules  by  HF  mol¬ 
ecules  gei\erally  brings  into  existerKe  the  following  kinds  of  defects: 

1 .  HF  itself  in  lattice  points.  We  denote  their  concentration  as  Hhel-  The  sequence  of  let¬ 
ters  HFL  implies  that  the  proton  (H)  and  L-defect  (L)  are  still  bound  with  fluorine. 

2.  HF  without  the  L-defech  «hf- 

3.  HF  without  a  proton;  npL- 

4.  HF  without  the  L-defect  and  a  proton:  rif. 

5.  H3Cri  ion;  Mh30^- 

6.  OH"  ion:  now" 

7.  L-defect; 

8.  D-defecb  «[> 

For  these  eight  unknown  variables  to  be  found,  eight  equations  are  required.  Three  of  those  are 
evident 


^HF  =  ”HFL  +  ”HF  ”FL  +  "F 


where  Nhf  denotes  total  concentration  of  impurity  molecules  HF. 

The  next  equation  implies  that  H30^  ions  are  formed  either  coupled  with  OH"  ions  from 
neutral  water  molecules,  or  when  a  proton  is  released  from  the  HF  molecule 


«H30+  =  «oh-  +  ”F+«FL-  (120) 

In  the  same  fashion,  L-defects  can  appear  either  coupled  with  EKlefects  from  a  pair  of  nor¬ 
mal  hydrogen  bonds,  or  can  be  released  from  the  HF  molecule.  Therefore 

«L  =  ”D  +  ”HF  +  ftp  •  (121) 


The  rest  of  the  five  equations  will  arise  from  the  law  of  mass  action  being  applied  to  revers¬ 
ible  reactions,  resulting  in  the  apprearance  of  charge  carriers.  So,  from  the  reaction 

HF  +H20«F-h-H30+-hL 

it  follows  that 

”H30^.'  =  Koi  «  expC-^Voi/fcaT)  (122) 

«HFL 

where  Kqi  is  the  dissociation  constant  of  this  reaction  and  Wqi  is  the  activation  energy.  Let  us 
write  down  equations  for  the  other  four  reactions 


(^Hao'*’  ■  ”fl)  /”hfl  =  “  exp("’^iF/*BT) 

(”HF  •  ”l)  /”HFL  =  f^OF  exp(~^OF/*^B^) 
(«F  •  ”l)  /”fl  =  f^i,o  “  exp("^i-o  /*bt) 


(123) 

(124) 

(125) 
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("F  •  «H3o")  /"HF  =  ^0,1  “  exp^-i^oj/tBr)  (126) 

If  the  dissociation  constants  are  known,  then  all  concentrations  can  be  found.  And,  if  activatim 
energies  are  known  in  addition,  then  the  temperature  dependence  of  all  concentrations  can  be 
predicted.  Of  course,  even  though  this  is  an  empirical  (rather  than  a  strictly  theoretical)  consid¬ 
eration,  the  actual  calculations  are  rather  cumbersome.  Sometimes  it  is  better  to  simplify  the 
consideration,  essentially  taking  into  account,  for  instance,  that  the  dissociation  energy  of  the 
L-defect  released  from  the  HF  molecule  is  extremely  small,  and  hence  fluorine  atoms  cannot 
keep  back  the  L-defects 


«F  =  NhF  +  «D  •  (127) 

If  in  addition  Nyjp  concentration  is  much  greater  than  the  intrinsic  equilibrium  concentration 
of  Bjerrum  defects,  then 


«L  =  Mhf- 


(128) 


In  this  particular  case,  finding  H30‘*^  ion  concentration  will  also  be  simplified,  since 

«HF  =  WhF  -  «F  =  «H30'^-  (129) 


Then,  from  eq  126 

«H30^  =  i[-Ko,I±(Ko,l"  +  4NHF-  (130) 

where  we  must  take  the  +  sign. 

At  high  temperatures  or  low  concentrations  of  HF 

(^o,i>>('(hf 

and 

”H30'^  =  N  hF  • 

At  low  temperatures  or  high  concentrations  of  HF 
0,1  <  <  ^  HF 

and 

”H30'^  =  [1^0,1  ■  (^hf]  '  ■ 

EXPERIMENTAL  TECHNIQUES  FOR  INVESTIGATION  OF 
CONDUCTIVITY  AND  DIELECTRIC  PERMITTIVITY  OF  ICE 

Measurements  of  the  electrical  properties  of  ice  can  provide  useful  information,  not  only  in 
respect  to  the  a  and  e  of  ice,  but  about  impurities  and  even  the  age  of  ice  as  well.  In  addition, 
since  the  same  defects  (protonic  charge  carriers)  are  responsible  for  both  electrical  and  mechani- 


(133) 

(134) 


(131) 

(132) 
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Osdllator 


Figure  18.  Schematic  sketch  ofac  bridge  circuit  to  measure  ac 
impedance  of  ice  sample,  Z  =  R+  lAtoC. 

cal  properties  (plasticity,  anelastic  relaxation  time),  electrical  measurements  in  ice  appear  to  be 
of  use  to  diverse  scientists  working  on  this  important  material.  Hierefbre,  it  is  necessary  to  get 
an  idea  of  basic  experimental  techniques  applied  in  investigations  of  a  and  e. 

Measuring  circuits 

The  most  common  measuring  circuit  for  electrical  measurements  in  ice  is  the  ac  bridge,  shown 
schematically  in  Figure  18.  Its  operational  principle  is  clear  from  the  sketch.  The  bridge  is  ad¬ 
justed  until  the  minimal  readings  of  the  null  indicator  are  observed.  The  ac  bridges  actually  used 
are  both  privately  made  and  industry  manufactured.  A  proper  selective  amplifier  or  even  a  lock- 
inamplifier  can  play  the  role  of  null  detector.  Such  circuits  are  known  to  be  used  in  file  frecjuency 
range  from  lO"^  up  to  10^  Hz.  However,  the  use  of  ac  bridges  at  frequencies /<  1  Hz  encoun¬ 
ters  ccHisiderable  difficulties.  As  a  rule,  this  circuit  is  employed  in  the  frequency  range  1-10^  Hz. 
With  the  aid  of  an  ac  bridge,  the  equivalent  electrical  resistance  and  capacitance  C(f)  (con¬ 
sidered  to  be  in  parallel)  of  ice  is  measiued  at  a  fixed  fi«<iuency.  Transformation  of  R(f)  and  C(f) 
into  a{f)  and  8(/)  is  not  difficult 

o(/)  =  ^;e(/)  =  ^C(/)  (135) 

sR\f) 
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Figure  19.  Schematic  representation 
of  "loop-method" for  ice  ac  impedance 
measurement. 


where  L  and  s  are  specimen  ttuckness  aivd  surface  area. 

Another  extensively  used  dicuit  is  represented  in  Hgure  19.  This  is  the  so-called  1(X^  meduxi 
An  XY  recorder  can  be  used  as  an  "osdUoscope"  in  the  low-frequency  range  (f^  1  Hz).  The 
circuit  has  proven  itself  to  good  advantage  in  the  frequency  range  of  1(H  to  10^  Hz. 

It  can  be  easily  shown  that 

R{f)  =  RoVo/VVa<2-Vi2;  C{f)  =  ^  (136) 

Vo  InfRo 

where  Vg  is  the  amplitude  of  the  oscillator  voltage.  From  here  a  and  e  can  be  evaluated  using 
eq  135.  It  is  necessaiy  for  measuring  resistance  that  Rq  be  much  smaller  than  ice  impedaitce  over 
all  frequency  ranges  used  in  the  measurements. 

The  third  technique  is  based  uptm  measurement  of  current  passing  throu^  the  sample  as  a 
function  of  time  I(t),  after  a  rectangular  voltage  has  been  applied.  So  that 

V(f)=j°'^^°  (137) 

|Vo  =  EL,tkO 

In  that  case,  whoi  the  electrodes  used  ensiue  an  efficient  charge  exchange  with  ice  (a  space 
charge  biiilt  up  in  ice  is  absent  near  the  electrodes),  a  dependence  I(t),  shown  in  F%ure  3,  is 
observed.  As  can  be  seen  from  dus  figure,  and  Og  are  determined  by  the  initial  and  final  (t 
»^d)  values  of  the  current  respectivdy.  Ihe  hig^-frequency  capacity  of  a  specimen  = s€(fi^/ 
L  defines  a  short  time  pulse  of  the  current  (not  shown  cm  Fig.  3)  with  characteristic  decay  time, 
T  =  Co,  R^  where  is  the  resistance  of  the  measuring  circuit. 


Static  dielectric  permittivity  can  easily  be  found  from  the  relationship 

os 

eoes  =  4-  f[l(t)-lH]dt.  (138) 

sV  J 
0 

All  three  methods  can  be  easily  computerized,  provided  that  digital  storage  oscilloscopes  and 
lock-in  amplifiers  operated  by  computers  are  used. 

Electrodes 

The  fabrication  of  proper  electrodes  for  ice  encounters  great  difficulties,  since  all  metals  are 
electronic  conductors,  while  ioe  is  a  protonic  axiductor.  In  many  cases  ffie  use  of  improper  elec¬ 
trodes,  causing  only  partial  charge  exchange  with  ice,  leads  to  observation  of  "false"  extra  dis- 
persiorrs  and  to  distortion  of  the  results  of  low-frequervry  measurements.  The  following  types 
of  electrodes  used  with  ke  can  be  foimd  in  the  literature: 

1.  Blocking  electrodes  (Gross  1975,  Mounier  and  Sixou  1969,  Zaretskii  1991). 

2.  Cold  plasma  (Auvert  and  Kahane  1973). 

3.  Palladiiun,  or  palladivun  black,  saturated  with  hydrogen  (BuUemer  et  al.  1%9,  Petradco 
etal.1983). 

4.  Ohmic  electrodes  with  double  electric  layers  at  the  ice/metal  interfece,  (Petrerrko  and 
Chesnakov  1990a,b/:). 

5.  Liquid  electrodes  (von  Hippel  et  al.  1973). 

6.  Ion  exchange  membranes  (Kahane  1969). 

Every  electrode  type  listed  above  has  been  used  in  many  experiments.  The  references  given 
indicate  either  the  first  lise  of  such  electrodes  or  the  most  detailed  description  of  the  techruque. 
All  electrodes  listed  yield  good  reproducible  results,  but  from  our  perspective,  the  most  reliable 
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and  easOy  built  are  blocking  electrodes,  in  whkh  the  charge  exchange  is  totally  eliminated  by 
thin  dielectric  films  placed  between  ice  and  metal,  atul  the  Ohmic  electrodes  described  in 
Petrenko  and  Chesnakov  (1990a,b,c),  in  which  the  charge  exchange  is  perfect,  siiKe  ekctrcmk 
energy  levels  in  ice  and  metal  become  equal  (see  the  Charge  Exchange  at  Ice/Metal  Interfaces  sec¬ 
tion).  Still,  it  is  probably  a  matter  of  taste. 

Surfoce  conductance  and  guard  rings 

As  we  will  see  later,  at  temperatures  T  >  -25°C,  even  in  pure  ice  a  considerable  surface  ccm- 
ductivity,  which  under  certain  conditions  exceeds  the  bulk  static  conductivity  of  ice,  is  observed. 
There  are  several  ways  of  checking  surface  ccmductance  (use  of  the  four-point  mehcxi,  for  in¬ 
stance  Qaccard  1966]).  However,  the  use  of  a  guard  ring,  sketched  in  Figure  20,  is  the  most  reli¬ 
able  and  widely  used  technique.  The  main  point  in  using  guard  rings  is  that  surface  currents 
will  not  enter  the  measuring  circuit,  thus  flowing  down  to  ground.  For  tinis  to  occur,  the  elec¬ 
tric  field  vector  must  be  parallel  to  the  surfoce  at  every  point  (to  prevent  surface  currents  from 
flowing  down  into  the  bulk  and  "back").  It  is  important,  hence,  that  the  potential  of  the  guard 
ring  be  almost  he  same  as  at  the  bottom  electrode.  The  last  condition  is  provided  when  R  (in 
Fig.  20)  is  very  small  compared  to  the  ac  impedance  of  the  ice  sample. 

Influence  of  inhomogeneity  on 

the  frequency  dependence  of  ice  conductivity 

As  a  final  remark  on  experimental  techiuques,  let  us  consider  the  influence  of  the 
inhomogeneous  distribution  of  charge  carriers  on  the  shape  of  the  frequency  dependence  of 
conductivity.  For  homogeneous  ice,  complex  conductivity  a  deperrds  on  circular  frequency  as 
(eq67) 


1  -i  ano 

where,  usually,  »  Og.  Using  eq  71 

5((o)  =  — (139) 
(l-r<inD)<t» 


If  a  depends  on  r  a  =  a(co,j:),  where  r  is  the  axis  perpendicular  to  he  electrodes,  hen  thin  lay¬ 
ers  wih  thickness  dx  will  be  cormected  in  series  and  he  "apparent"  conductivity  ofap  will  be 


1 


'  -1 
dx 

o(oj,  x) 


<t»  (l-toTiW 


(140) 


where 


L 

^DM  =  f  (141) 

^  0 

HerKe,  we  have  again  in  eq  140  a  regular  Debye's  dispersion,  as  in  eq  139,  but  wih  an  "aver- 
age"  Tq.  On  he  oher  hand,  when  he  inhomogerveity  of  5(0))  exists  along  he  electrodes 
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®ap  =  g  2)  dy  dx  (142) 

S 

and  Cap  is  the  sum  of  dispersions  with  different  Tq  (and  (Od)  and  of  different  amplitudes  (o«  - 
Og).  Indeed,  in  a  real  experimental  situation,  ice  samples  prepared  for  electrical  measurements 
very  often  have  an  inhomogeneous  distribution  of  charge  carriers.  That  ii\homogeneity  may  arise 
from  the  presence  of  impurity  concentration  gradients,  dislocations  arwl  so  on. 


REVIEW  OF  EXPERIMENTAL  RESULTS  ON 

ICE  CONDUCTIVITY  AND  DIELECTRIC  PERMITTIVITY 

The  very  first  scientific  research  in  ice  COTiductivity  of  which  this  author  is  aware  began  more 
than  100  years  ago  (Ayrton  and  Perry  1877).  At  present,  the  list  of  papers  devoted  to  the  electri¬ 
cal  properties  of  ice  contains  several  hundred  items.  It  is  extremely  difficult  to  conduct  a  com¬ 
parative  and  critical  review  of  such  a  vast  amoimt  of  material.  It  was  already  a  huge  task  in  the 
early  seventies,  when  Hobbs  (1974)  was  working  on  his  famous  book.  Since  then  the  amount 
of  information  in  this  field  has  more  than  doubled,  and  here  we  present  the  material  in  a  com¬ 
pletely  different  way. 

First  of  all,  notice  that  all  theoretical  and  model  considerations  introduced  previously  in  this 
report  are  based  upon  critical  analysis  of  the  totality  of  experimental  data.  Therefore,  our  pri¬ 
mary  goal  now  is  to  illustrate  the  validity  of  these  ideas  in  the  context  of  some  particular  ex¬ 
perimental  research  and  to  provide  the  reader  with  some  useful  numerical  data  on  the  abso¬ 
lute  values  of  conductivity,  activation  energies,  charge  carrier  mobility,  relaxation  times  and 
dielectric  permittivity,  and  the  values  of  electric  charges. 

In  pursuing  this,  inasmuch  as  the  amount  of  material  available  is  much  more  than  is  neces¬ 
sary  for  illustration,  in  the  future  we  will  be  guided  by  the  following  criteria.  We  will  consider 
first  mciinly  recent  works  not  included  in  Hobbs'  book  and,  secondly,  works  that  are  remark- 


Figure  21.  Frequency  dependences  of  pure 
ice  dielectric  permittivity  e  and  conductiv¬ 
ity  o;T=  270.2  K.  The  letter  "D"  indicates 
Debye  dispersion  (after  Noll  1978). 
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able  in  some  asp>ect  (original  techniques,  precision  of  results,  tfie  role  they  played  in  establish¬ 
ing  modem  understaiuling  of  the  subject,  etc.).  Of  course,  the  choice  is  still  rather  arbitrary,  and 
the  author  presents  his  apologies  in  advance. 

We  will  begin  our  review  by  discussing  the  general  form  of  dispersion  curves  a(<B)  and  e(o)). 
As  we  have  seen  earlier,  Jaccard's  theory  predicts  the  presence  of  only  one  singularity  in  the 
spectrum  (single  dispersion),  as  shown  in  Figures  2  and  5.  Such  dispersions  cure  observed  when 
a  limited  rang^  of  hequendes  close  to  tOo  is  used.  When  dre  frequency  range  is  widerred,  an  extra 
dispersion  appears.  All  four  low-frequency  dispersions  (including  Debye  dispersion)  predicted 
by  Jaccard's  model  and  modified  for  a  finite-dimension  specimen  were  detected  by  Maidique 
et  al.  (1970)  in  ice  specimens  being  annealed  for  10  months. 

Three  separate  dispersions  for  co  <  (Oq  were  reserved  by  Noll  (1978)  (see  Fig.  21).  Many  other 
authors  reported  the  observation  of  one  extra  dispersion  at  fiequordes  less  than  <i)j>  As  we  have 
seen  in  the  Electrical  Properties  of  Ice  of  Finite  Size  sectiorv  these  disp>ersions  are  caused  by  screen¬ 
ing  layers  formed  near  the  ice  specimen  siuface  that  consist  of  various  charge  carriers.  As  men¬ 
tioned  already  in  that  section,  such  dispersions  can  yield  a  lot  of  mforrruition  about  concentra¬ 
tion  and  mobility  of  charge  carriers  (Zaretskii  et  al.  1987a,  2^aretskii  1991).  Since  these  disper¬ 
sions  appear  at  the  "electric"  boundaries  of  the  ice  specimen,  where  an  electric  potential  drop 
occurs,  then,  if  electrodes  providing  a  perfect  charge  exchange  with  ice  are  used,  these 
singularities  disappear  at  frequencies  oa  <  ojd- 

Many  authors  also  note  the  presence  of  small  additional  dispersions  in  some  ice  specimens 
at  frequencies  <n  >  oid-  Maidique  et  al.  (1970),  for  instance,  reported  observing  two  such  disper¬ 
sions.  The  nature  of  these  high-frequency  dispersions  is  still  imclear  and  under  discussion.  It  is 
possible  that  all  (or  at  least  part)  of  these  dispersions  are  determined  by  the  existence  of  local 
inhomogeneties  in  ice,  which  we  discussed  in  the  previous  section. 

Measuring  dispersion  curves  a(a))  allows  us  to  determine  the  values  of  static  (Oj)  and  high- 
frequerKy  (o J  conductivities.  From  measuring  o„  and  Oj  at  various  temperatures,  the  activa¬ 
tion  energies  Eg*  and  Egg  c^n  be  calculated.  Since  a  «  «  •  |i,  these  activation  energies  come  out 
to  be  the  sum  of  1)  activation  energies  of  the  creation  of  charge  carriers  (E^b  and  Eaj  for  Bjerrum 
defects  and  ions  respectively)  and  2)  activation  energies  of  defect  motion  E^j-  Experiments  car¬ 
ried  out  on  ice  doped  with  various  impurities  allow  us  to  predict  charge  carrier  concentrations 
and  thus  to  separate  values  of  n  and  p  and  to  separate  the  activation  energies  of  creation  and 
motion  of  charge  carriers  as  well.  The  results  of  such  measurements  are  summarized  in  Tables 
1-3.  We  present  below  the  "guidance"  for  these  experimental  results. 

nie  most  well  established  conductivity  parameters  for  piue  monocrystalline  ice  are  evidently 
high-frequency  conductivity  cs„  and  its  activation  energy  Ea„.  Good  reproducibility  of  these 
quantities,  determined  by  the  conductivity  of  m^ority  charge  carriers  (L-defects  for  T  >  -50°C), 
is  attributable,  first,  to  high-frequency  measurements  not  depending  upon  electrode  type.  Sec¬ 
ond,  it  appears  to  be  comparatively  easy  to  attain  a  degree  of  purity  in  ice,  such  that  is  de¬ 
termined  by  intrinsic  carriers  and  does  not  depend  upon  residual  impurities.  The  same  holds 
for  Tq  and  E.^  since,  according  to  eq  69  and  70,  Xp  «  always  when  »  Og,  i.e.,  when  the 
majority  type  of  charge  Ccirrier  is  clearly  distinguished.  These  four  values  at  T  =  -10°C  are  as 
follows 


Ea~  =  (o.58  ±  0.03)  eV;  Ex  =  (0.58  ±  O.Ol)  eV 
=  (l.8  ±  0.2)  10"®  a~^  m-l ;  td  s  5.10“®s 


(143) 


The  Debye  dispersion  is  reproduced  very  well  in  perfect  crystals — pure,  homogeneous  and 
annealed  (for  relaxation  of  internal  stresses). 
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Table  1.  Data  on  high-frequency  (o^  and  static  (o,)  conductivities  of  ice. 
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KOH  doping  -23  to -203  -2xl0-‘  10^  0  0.4  Howe  and 

0.032  mol/L _ (-23°C)  (-23°C)  g>-100°C)  (T>-100°C) _ 0 _ Whitworth  (1989) 


Considerably  greater  scattering  is  observed  in  data  from  measurements  of  static  conductiv¬ 
ity  Os  and  its  activation  energy  E^.  Values  of  appear  to  lie  between  0  and  0.7  eV  and  Oj  (- 
lO'C)  varies  in  limits  from  10^  down  to  6  x  10“'*^  m“' .  This  is  explained  by  the  difficulty  in 

purifying  ice  to  a  level  at  which  intrinsic  ions  dominate  over  ions  supplied  by  impurities.  That 
is  why  von  Hippel  (1971)  and  von  Hippel  et  al.  (1973)  drew  the  conclusion  that  there  are  no  in¬ 
trinsic  ions  in  ice  and  that  all  ions  originate  in  impurities.  If  we  refer  to  Table  I  from  their  work 
(von  Hippel  et  al.  1971)  for  the  value  of  Oj  (2.5  x  10^  £1“^  m“' )  and  compare  it  with  the  minimal 
observed  value  6.4  x  10“^®  then  we  might  agree  widi  them — their  ions  were  of  impurity  ori¬ 
gin.  The  question,  however,  remams  open:  Were  the  specimens  in  the  experiments  of  Worz  and 
Cole  (1969)  or  Petrenko  et  al.  (1983)  pure  enough  that  their  conductivity  may  be  considered 
intrinsic?  The  following  general  trend  is  observed:  The  more  pure  the  ice  and  the  smaUer 
the  greater  is  E^^  Therefore,  the  value  of  the  activation  energy  of  creation  of  an  ion  pair  E^  - 
0.96  eV  (Eigen  and  De  Mayer  1956, 1957, 1958;  Eigen  et  al.  1964),  considered  for  a  long  time  to 
be  commorvly  acknowledged  truth,  should  be  treated  as  the  lowest  possible  limit  of  this  quan¬ 
tity.  As  the  "available"  upper  limit  we  can  take  £;„  =  2£a5  =  1 .4  e V  from  the  work  of  Petrenko  et 
al.  (1983). 

Good  experimental  confinnation  of  theoretical  concepts  about  the  effect  of  impurities  on  ice 
conductivity  (the  Protonv:  Charge  Carriers  Introduced  by  Doping  section)  has  been  indicated  in 
many  papers  (see,  for  instance,  Takei  and  Maeno  1984, 1987;  Camplin  et  al.  1978).  The  most  ef¬ 
ficient  impurities  for  the  increase  of  H30‘''  ion  and  L-defect  concentrations  appear  to  be  the 
adds  HE  and  HCl.  KOH  doping  increases  OH“  ion  and  L-defect  concentrations,  so  diat,  firstly, 
ice  with  amazingly  large  conductivity  (up  to  =  10  m~^  at  T  =  CPC)  can  be  obtained  and, 

secondly,  considerable  conductivity  (a^)  is  observed  at  temperatures  as  low  as  -200°C  (Zaretskii 
etal.  1988). 

The  determination  of  charge  carrier  mobilities  in  ice  is  an  extremely  difficult  task,  since,  be¬ 
cause  of  small  values  of  pj  and  the  "hop"  mechanism  of  motion,  the  regular  technique  for  mea¬ 
suring  the  HaU  effect  is  inapplicable  in  this  case.  (The  Hall  effect  will  be  discussed  in  the  report. 
Surface  of  Ice,  listed  in  the  Foreword.)  Hopes  of  using  "saturation  currents"  for  determining  pj  are 
not  justified,  because  of  reasons  pxrinted  out  by  Maidique  et  al.  (1971)  in  their  critique  of  Eigen's 
work,  and  also  because  when  high  voltage  is  applied  to  a  system  with  ice/ metal  interfaces, 
various  physical  phenomena  that  make  the  interpretation  of  results  more  complicated  are  ob¬ 
served  (see  Part  II). 

Nonetheless,  several  new  techniques  for  determining  charge  carrier  mobilities  in  ice  have  been 
developed  and  employed  during  the  last  25  years.  These  are  the  stationary  and  transient  injec¬ 
tion  currents  measurement  technique  (Eckener  et  al.  1973,  Petrenko  et  al.  1983,  Petrenko  and 
Ryzhkin  198^),  the  ice  field  effect  transistor  (Petrenko  and  Maeno  1987),  and  the  low-fi-equency 
dispersion  measurement  technique  (Petrenko  and  Ryzhkin  1984a,  Zaretskii  et  al.  1987,  Zaretskii 
1991).  The  idea  of  determining  pj  and  from  the  parameters  of  low-frequency  dispersions  was 
presented  in  the  Electrical  Properties  of  Ice  of  Finite  Size  section.  Other  techniques  mentioned  above 
will  be  described  and  discussed  later  in  Part  II.  The  results  of  determining  charge  carrier  mo¬ 
bilities  in  ice  are  summarized  in  Table  2.  At  first  glance  only  the  data  on  L-defect  mobilities  seem 
to  be  satisfactory.  At  T = -10°C,  P4 = (2-5)  x  ICf^  m^/ V  s.  This  mobility  drops  exponentially  with 
decreasing  temperature,  the  activation  energy  being  0.19-0.235  eV  (see  Table  1). 

StiU,  there  is  no  reliable  mformation  about  D-defect  mobility.  There  is  only  indirect  evidence 
thatp3<p4. 

Data  on  icm  mobilities  are  distinguished  at  first  glance  by  large  scattering  (from  5  x  10"^  down 
to  25  x  10^  for  H30‘*'  and  from  4  x  10“^  down  to  2  x  10”^°  for  OH“).  However,  if  we  discard 
large  values  obtained  by  the  now-discredited  method  of  saturation  currents  and  values  obtained 
through  measurements  of  low-ltequency  dispersions  (which  still  are  not  accuracte  enough)  as 
well,  it  appears  that  positive  ion  mobility  lies  within  the  limits  2.7  x  10^  and  1.1  x  10“^  m^/V  s 


39 


Table  2.  Mobilities  of  charge  cairieis  in  ice. 


T 

(°C> 

Ml 

(m^/Vs) 

M2 

(m^/V  s) 

Ml 

(m^/Vs) 

Mi 

(rn^/Vs) 

Method  used 

Reference 

-13  to -36 

(l.l±0.1)xl0-^ 

Analysis  of 

conductivity 

decay 

Kunst  and 
Warman  (1983) 

-10 

2.5xl0-» 

o,and 

measurements 

Wdrz  and  Cole 
(1969) 

Oto-20 

5x10-^ 

Saturation 

currents 

Maidique  et 
al.  (1988) 

-123 

S4xl0-^ 

Doping  with 

KOH 

Zaretskii  et 
al.  (1988) 

-5to-40 

i(r®- 10-^ 

Proton  injection 

PetrerJto  et  al. 
(1983) 

0 

2.7xl0-» 

5xl0-» 

Doping  with 

HP 

Camplin  et  al. 
(1978) 

-10 

1.8x10-* 

Doping  with 

HP 

Jaccard  (1959) 

:^xio-“> 

Howe  and 
Whitworth 
(1989) 

-10 

(2.4±1.6)xir^ 

Saturation 

currents 

Bullemer  et  al. 
(1969) 

-10 

3.5x10-^ 

HP  doping 

Camplin  and 
Glen  (1973) 

^.1 

9.2xl0-» 

2.7x10-* 

Ice  field 
transistor 

Petrenko  and 
Maeno  (1987) 

-14510-178 

(9±l)xlO^ 

Transient 
injection  of 
protons 

Eckener  et  al. 
(1973) 

-10 

=10-* 

1.7x10-* 

Low-frequency 

dispersions 

Zaretskii  et 
al.  (1987a) 

-33 

(6.7±0.8)xl0-^ 

Low-frequency 

dispersions 

Zaretskii  (1991) 

-20 

6x10-* 

Proton  injection 

Petrenko  and 
Ryzhkin  (1984b) 

-10 

2.4x10-* 

Recombination 

ir^tion 

Petrenko  and 

Chesnakov 

(1990c) 

— not  bad  at  all!  The  author  of  this  review  prefers  the  data  obtained  in  transient  injection  ex¬ 
periments  (Eckener  et  al.  1973),  in  whidx  a  minimum  of  assumptions  was  used.  The  same  H3O* 
ion  mobility  was  found  by  Petrenko  aivd  Maeno  (1987)  using  field  effect  transistors.  So,  it  is  likely 
that  Pi  =  9  X  1(1^  m^/ V  s.  We  can  consider  it  Scifely  established  that  botti  H3Cy'  and  OH"  ions 
move  with  zero  activation  energy;  this  is  in  good  agreement  with  concepts  about  the  tunneling 
of  protons  along  hydrogen  bonds. 
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Table  3.  Ratio  of  effective  electric  charge  of  H3O'*’  ion  (ci)  and 
of  D-defect  iej}  to  proton  charge  (e). 

ei/e  eye  Reference 


(0.61  ±  0.02)  0.39  ±  0.01  Zaxetskii  et  aJ.  (1968) 

0.73  0.44  Camplin  et  al.  (1978)* 

0.62  ±  0.01  0.38  ±  0.01  Hubmann  (1979a) 

0.626  0.374  Take!  and  Maeno  (1987) 

0.64  ±  0.03 _ 0.36  ±  0.03 _ Scheiner  and  Nagle  (1983) 

*  Data  fitted  without  requiring  e]+e3  =  e;  may  cause  some 
errors  in  their  values  in  Tables  1-3. 


In  this  section  we  are  not  going  to  c(»i- 
sider  in  detail  the  surface  cc«\ductivity  of 
ice,  since  its  nature  evidently  differs  appre¬ 
ciably  from  Jaccard's  model,  which  is  actu¬ 
ally  a  dieoretical  framework  for  this  report 
It  will  be  more  convenient  to  discuss  the 
problem  of  the  surface  conductivity  of  ke  in 
one  of  our  future  reports,  specially  devoted 
to  ice  surfaces.  Here  we  restrict  ourselves 
to  Figure  22,  adopted  from  the  work  of 
BuUemer  and  Riehl  (1966),  which  illustrates 
die  relationship  between  bulk  static  conduc¬ 
tance  and  surfooe  conductivity. 

Let  us  now  briefly  consider  experimen¬ 
tal  results  on  measurements  of  ice  static  di¬ 
electric  permittivity  Eg  and  compare  these  results  with  the  predictions  of  Jaccard's  model. 

Note  first  of  all  that  theory  yields  quite  a  consistent  dependeiKe  of  e  as  a  function  of  (o,  both 
in  the  region  of  principal  Debye  dispersion  (see  eq  74)  and  in  the  low-frequency  dispersion  re¬ 
gion,  as  was  described  earlier.  The  predicted  absolute  value  of  static  dielectric  permittivity  = 
(-10®C)  and  its  increase  with  temperature  decrease  (see  eq  72)  are  also  confirmed.  As  we 
wiU  see  in  Part  n,  the  value  of  e,  is  brilliantly  predicted  in  those  cases  when  the  values  of  Oi  and 
Ob  become  comparable  (so  called  "crossover,"  discussed  in  Part  II),  which  again  follows  from 
eq  72.  Yet,  the  theory  in  its  present  form  has  some  problems.  For  example,  it  does  not  account 
for  the  sligjit  anisotropy  of  e  observed  in  a  number  of  experiments.  Let  us  refer  to  one  of  those 
recent  works  by  Takei  and  Maeno  (1987)  where  they  found  that 


Figure  22.  Temperature  dependence  of  bulk  and 
surface  conductance  cf  pure  ice  specimen  (after 
Bullemerand  Riehl  1966). 


es-  e- 


2.41  X  lO^K 
T-2.2K 


(144) 


Es- 


2.3  X  lO^X 
T-  48.7 K  ' 


(145) 


As  can  be  seen  Iran  eq  144  and  145,  instead  of  the  predicted  behavior  of  (Cg  -  O  according  to 
Curie's  law  (see  eq  72) 

(%-£<»)“- (146) 
♦  i 
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the  Curie-Weiss  law  is  observed 


(£s 


_J _ 

T-Tc 


(147) 


which  implies  a  transition  to  the  ordered  phase  at  T  S  T^.  There  is  a  coi\siderable  discrepancy 
in  Tj.  in  the  data  of  various  authors  who  performed  experiments  on  pure  mtwrocrystalline  ice. 
It  ranges  from  0  K,  according  to  Ruepp  and  (1%9),  up  to  48.7  K,  according  to  Takei  and 
Maeno  (1987).  Some  errors  in  measuring  in  "pure"  ice  at  low  temperatures  can  arise  because 
of  increasing  screening  lengths  and  an  apparent  drop  of  (see  the  Electrical  Properties  of  Ice  of 
Finite  Size  section).  Hence,  critical  temperatures  observed  in  pure  ice  can  be  treated  as  the  lower 
limit  of  actual  critical  temperatures.  In  fact,  Tajima  et  al.  (1982)  fouird  the  transition  of  KOhi- 
doped  ice  into  a  partially  ordered  state  to  occur  at  T  =  72  K . 

In  general,  we  should  infer  that  at  present  Jaccard's  theory  (along  with  its  further  develop¬ 
ment)  describes  well,  both  quantitatively  and  qualitatively,  the  diverse  electrical  properties  of 
ice.  As  a  brilliant  illustration  of  this,  1  ask  the  reader  to  admire  Figure  21  once  again! 
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PART  II.  ADVANCED  TOPICS  AND  NEW  PHYSICAL  PHENOMENA 


In  this  part  we  treat  more  advanced  and  special  questions  concerning  the  electrical  proper¬ 
ties  of  ice.  Most  of  these  belong  to  the  frontiers  of  research  in  the  "physics"  of  ice,  as  an  object  of 
solid  state  physics.  The  results  and  physical  models  presented  here  should  be,  for  the  most  part, 
of  interest  to  physicists  and  chemists  doing  direct  research  on  ice,  rather  than  the  larger  group 
of  readers  who  intend  to  use  informahon  about  ice  properties  only  for  some  applications. 


CHARGE  EXCHANGE  AT  ICE/METAL  INTERFACES 

The  electric  charge  exchange  between  ice  and  metal  determines  the  conditions  for  electric 
current  flow  through  ice  via  metallic  electrodes.  Such  an  exchange  plays  an  essential  role  in  the 
process  of  chemical  reactions  at  the  ice/metai  interface.  Common  metal  oxidation  can  serve  as 
an  example  of  such  reactions.  A  future  report  devoted  to  electro-optical  effects  in  ice  finds  that 
photochemical  reactions  also  proceed  at  such  interfaces,  which  implies  electron  exchange  be¬ 
tween  ice  and  electronic  conductors  (metals,  semiconductors). 

It  has  been  common  knowledge  for  a  long  time  that  stable  electric  currents  cannot  pass 
through  an  ice/metal  interface,  even  when  precious  metals  (Pt,  Au),  which  have  a  surface  that 
is  not  oxidized  by  contact  with  ice,  are  used.  Applying  constant  voltage  to  ice  with  two  metal¬ 
lic  electrodes  results  in  a  short  current  "pulse,"  which  attenuates  exponentially  according  to  the 
presence  of  several  relaxation  times  (see  the  Electrical  Properties  of  Ice  of  Finite  Size  section  and 
the  previous  section).  That  is,  no  regular  metallic  electrodes  are  capable  of  providing  a  stable 
charge  exchange  between  ice  and  metal,  which  is  necessary  for  a  stable  dc  flow. 

Hydrogen  saturated  Pd  electrodes  are  an  exception,  but  in  this  case  Pd  and  ice  exchange 
protons,  rather  than  electrons  (see  the  Proton  Injection  from  Pd  Electrodes  into  Ice  section  for  de¬ 
tails).  Note  also  that  there  is  a  charge  exchange  between  metal  and  a  quasi-liquid  surface  layer, 
since  there  are  surface  currents  present. 

There  are  no  obstacles  to  dc  flow  through  water  and  two  metallic  electrodes.  Insofar  as  elec¬ 
tronic  level  structures  in  water  and  ice  are  rather  similar  (which  can  be  inferred  from  optical 
measurements),  such  differences  in  the  properties  of  ice/metal  and  water/metal  interfaces  are 
puzzling. 

The  key  to  the  mystery  was  foxmd  at  first  in  the  work  of  Evtushenko  et  al.  (1988),  in  which 
for  the  first  time  a  technique  was  developed  that  ensures  a  sfable  dc  flow  through  the  ice/metal 
interface.  The  ultimate  understanding  was  achieved  thanks  to  work  by  Petrenko  and  Chesnakov 
(1990b)  and  Chesnakov  (1990).  I  present  below  a  brief  accovmt  of  the  results  of  these  investiga¬ 
tions. 

Since  the  electronic  level  structure  of  water  and  ice  are  similar,  we  will  start  with  water,  whose 
charge  exchange  with  metal  is  well  known  (see  for  instance  Bockris  1979,  Vihj  1973).  The  mu¬ 
tual  arrangement  of  electroruc  energy  levels  in  metal  electrodes  and  water  before  {V  =  0)  and 
after  {V  >  123  V)  a  voltage  was  applied  is  shown  in  Figure  23.  The  charge  exchange  being  real¬ 
ized  is  between  the  Fermi-level  of  the  metal  £f  and  the  hydrogen  level  £^2  of  water  at  the  cath¬ 
ode  (according  to  eq  23)  and  between  Ef  and  oxygen  level  E^  at  the  anode  (according  to  eq  24). 
The  "hydrogen"  system  includes  H30'^  ions,  and  atoms  and  molecules  of  hydrogen.  £1^2  can 
be  looked  upon  as  an  electron  chemical  potential  in  this  system.  The  same  is  true  for  ^02  and 
the  "oxygen"  system,  which  comprises  OH~  ions,  and  atoms  and  molecules  of  oxygen.  In  wa¬ 
ter  £h2  -  £02  =  A  =  123  eV.  For  efficient  electron  exchange  between  two  levels  to  be  realized, 
their  equalization  is  required.  In  the  case  of  water,  where  the  potential  difference  being  applied 
between  anode  and  cathode  V  >  1 .23  eV  =  A,  this  equalization  occurs  automatically  because  of 
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Figure  23.  Electronic  energy  levels  in  water  and 
two  metal  electrodes. 


the  double  diarged  layers  being  formect  it  v-o 

the  water/ metal  interface. 

The  double  layer  of  the  anode  consists 
of  positively  charged  holes  in  metal  and 
(negatively  charged)  OH~  ions,  "stuck"  to 
the  electrode  surface  in  water.  The  cathode 
double  layer  is  formed  by  electrons  in 
metal  and  ions  "stuck"  to  the  cath¬ 
ode.  The  action  of  such  charged  double  lay¬ 
ers  on  the  potential  is  analogous  to  the  ac¬ 
tion  of  flat,  charged  electrical  capacitors — 
they  provide  an  electrical  potential  drop, 
necessary  for  level  equalization.  If  V  >  A, 
the  rest  of  the  voltage  drops  at  water  bulk. 

Electrons,  passing  from  water  to  metal  and 
back,  overcome  an  energy  barrier  having  a 
width  about  that  of  the  double  charged 
layer.  In  the  case  of  water,  the  width  of  this 
barrier  is  several  angstroms  and  electrons 
penetrate  it  easily  by  tunneling.  In  the  case 
of  ice,  the  width  of  a  similar  layer,  which  is 
nothing  more  than  a  screening  layer  con¬ 
sisting  of  ions,  is  about  1  pm  (|im  =  lO^A), 
since,  first,  the  ionic  concentration  drops  by 

several  orders  of  magrutude  and,  second,  the  effective  dielectric  permittivity,  acting  at  a  range 
of  several  intermolecular  distances,  decreases  (from  =  80  down  to  3.2).  And  as  we  remember 
from  an  earlier  section,  the  screening  length  1/k  is  proportional  to  Ve^i  (see  eq  95,  we  have  to 
use  for  water  e  =  80  and  for  ice  e  =  e„  =  32).  As  a  result,  for  the  potential  difference  being  ap¬ 
plied  to  ice  with  two  metallic  electrodes  attached,  although  the  double  charge  layers  equaliz¬ 
ing  electronic  energy  levels  do  form,  the  width  of  these  layers  is  so  great  that  electrorts  carmot 
overcome  it  through  the  tunneling  process 
As  a  result,  dc  does  not  flow!  (A  =  123  +  0.06 
=  1.29  eV  is  also  greater  in  the  case  of  ice.) 

Here,  we  have  taken  into  account  the  latent 
heat  of  ice  melting  (0.06  eV). 

However,  as  was  shown  by  Evtushenko 
et  al.  (1988),  the  problem  of  charge  exchange 
between  ice  and  metal  can  be  solved,  pro¬ 
vided  that  at  first  thin  and  dense  charged 
double  layers  are  formed  at  the  metal/wa- 
ter  interface,  and  are  then  frozen  down  to 
T  <  -10°C.  These  thin  ion  layers  at  the  an¬ 
ode  and  cathode  provide  good  conditions 
for  charge  exchange  between  ice  and  metal 
electrodes.  The  layers  themselves  appear  to 
be  very  stable  and  can  exist  for  one  day  at 
T  =  -10°C  and  for  one  week  at  T<  -20'’C. 

Such  a  long  lifetime  is  explained  by  a  large 
discharge  time  and  constant  RC  of  such  a 
system,  where  C~10^  F/cm^  and  R  ^  10"^ 


U(V) 

Figure  24.  Current-^ltage  characteristic  of  ice  speci¬ 
men  with  two  Pt  electrodes.  Electrically  charge  double 
layers  were  formed  at  Pt/water  inteifaces  and  then  fro¬ 
zen  dawn;  T  =  -10°C  (after  Evtushenko  et  al.  1988). 
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Q  cm^,  which  yields  t  >  10*  s  (T  =  -10°C).  An  interesting  peculiarity  of  electrodes  processed  in 
this  way  is  that  they  let  current  pass  in  only  one  direction;  i.e.,  they  appear  to  be  quite  good  rec¬ 
tifiers,  as  indicated  in  Figure  24  (see  also  Fig.  30)  The  rectification  ratio  reaches  values  of  10*  and 
more.  Further  on,  we  will  refer  to  these  electrodes  as  ohmic,  because  of  their  abiUty  to  let  elec¬ 
tric  current  pass  essentially  without  voltage  drop  (above  =  13  V,  necessary  to  compensate  for 
"frozen"  potential  difference).  After  the  preparation  process  is  over,  the  electrodes  (with  the 
double  layers  of  ions  and  thin  layer  of  ice)  can  be  removed  from  the  bulk  of  ice  and  used  with 
other  specimens  (monocrystals  for  instance).  Utey  have  already  found  a  use  in  several  research 
works. 


RELAXATION  TIMES  OF  ELECTRIC  POLARIZATION 
AND  ELECTRIC  HELDS  IN  ICE 

The  question  of  electric  field  and  polarization  relaxation  rates  has  general  scope  and  can  be 
found  in  many  real  applications.  An  electric  polarization  can  arise  owing  to  both  various  exter¬ 
nal  effects  (external  electric  fields,  mechanical  loadings,  temperature  gradients)  and  internal 
changes,  such  as  crack  growth.  Once_electric  polarization  P  appears,  how  fast  (or  for  what  char¬ 
acteristic  time)  will  it  decay?  Since  P  will  be  "dissolved"  by  means  of  protonic  charge  carrier 
motion,  then  it  might  seem  we  have  already  found  an  answer  to  this  question  in  earlier  sections. 
First  of  all,  this  is  Debye  relaxation  time  Xq,  determined  by  the  majority  charge  carrier  motion. 
In  the  case  of  ice  having  boimdaries,  slower  relaxation  processes,  determined  in  terms  of  mi¬ 
nority  charge  carrier  motion,  can  be  observed  at  time  intervals  longer  than  Xp.  However,  we 
might  expect  that  the  better  part  of  the  polarization  will  disappear  with  the  most  rapid  charac¬ 
teristic  time  Xp.  Although  we  hear  this  answer  from  most  people  doing  electric  measurements 
on  ice,  it  is  absolutely  wrong! 

Let  us  try  to  illustrate  this  fac^Consider  a  very  large  ("infinite")  ice  volume,  and  a  relaxation 
process  of  electric  polarization  arising  in  a  finite  domain  of  this  ice  volume.  Obviously,  the 
following  equation,  familicir  from  electrodynamics,  holds 

div  £  =  — —  =  -div  P/eEq  (148) 

where  p  =  space  charge  density. 

Taking  into  account  that  curl  £  =  0  (quasi-stationary  approximation),  we  can  conclude  from 
eql48 

E=-£-  +  /{t)  (149) 

EEq 

where  the  vector  /  can  depend  only  upon  time.  But  taking  into  account  that  there  is  no  electric 
field  at  infinity  (that  is  why  we  consider  "infinite"  ice)  /  (f)  =  0  and 

£=-P/eeo.  (150) 

We  are  interested  in  this  particular  case  in  the  most  rapid  process,  determined  in  terms  of  ma¬ 
jority  charge  carriers,  i.e.,  carriers  having  maximal  conductivity.  In  the  case  of  ice  those  are  usu¬ 
ally  04  (L-defects) 

4  _ 

p  =1.^)/)=  -^3  /4  (151) 

/=! 
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(152) 


74=(-e3E 

t  1 

jj4{t)dt  =  74  .  (153) 

0 

Substituting  eq  150,  ^1  and  153  into  152,  we  obtain  the  desired  kinetic  equation,  defining  die 
time  dependence  of  £ 

E+Eh  =  0  .  (154) 

From  here 

E  =  E  (o)  exp  {-t  a)  (155) 

where  x  is  defined  as 


1  = 

(156) 

T 

tM  tD 

JL 

_  04 

(157) 

tM 

04 

(158) 

tD 

w(es-e«)' 

Note  that  Tq  is  already  familiar  to  us  as  Debye 
relaxation  time  (see  eq  43  and  71)  for  our  case  of 
CKie  type  carriers'  motion.  Xg^^  is  the  still  more  well- 
known  Maxwell  dielectric  relaxation  time,  intrin- 
sic  to  space  charge  "relaxation"  in  ordinary 
conductors.  As  can  be  seen  from  eq  156,  relaxation 
rate  (1  /x)  is  equal  to  the  sum  of  Debye  and  Max¬ 
well  relaxation  rates.  This  result  has  a  simple  in¬ 
terpretation.  The  driving  force  acting  upon  carri¬ 
ers  is  the  electric  field  in  the  Cc..e  of  Maxwell  re¬ 
laxation,  while  Etebye  relaxation  time  x^  is  asso¬ 
ciated  with  the  effect  of  the  configurational  vec¬ 
tor  Q .  In  gener^,  charge  carriers  in  ice  are  sub¬ 
jected  to  both  E  and  £2,  which  accelerate  re¬ 
laxation  processes. 

Of  course,  to  obtain  results  (eq  154-156)  it  is  not 
absolutely  necessary  that  the  ice  volume  be  "infi¬ 
nite."  Figure  25  illustrates  this.  For  ice  specimens 
cut  in  the  form  of  thin  plates,  so  that  P  initially 
was  normal  to  the  plates,  we  can  then  neglect  the 
end  effects  and  the  exact  solution  is  again  given 
by  the  equation 


Figure  25.  Characteristic  time  rwith  which 
electric  polarization  of  ice  relaxes  for  different 
boundary  conditions. 
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(159) 


P  +  P  /T  =0 

where  T  depaids  upon  forces  acting  upon  carriers,  which  depend  on  bourulary  ccnditions.  Note 
that  given  R^ce  as  an  ohmic  ice  resistance,  then  in  the  Figure  25c  circuit  we  have  for  Rjce/f^  -*  0 


H 

X 

X 


^  TM 


(160) 


and  forP,ce/R -»«» 


These  remarks  illustrate  the  importance  of  ice  resistance  and  input  circuit  resistance  ratio  for 
correct  determination  of  relaxation  times! 


RECOMBINATION  INJECTION  OF  IONS  INTO  ICE 

In  the  Concentration  of  Charge  Carriers  section,  we  considered  physical  laws  that  determine 
charge  carrier  concentration  in  ice  in  a  state  of  thermal  equilibrium  and  in  the  absence  of  exter¬ 
nal  excitations.  The  next  four  paragraphs  will  be  concerned  with  the  electrical  properties  of  ice 
in  strong  electric  Helds.  As  we  shall  see  below,  the  application  of  an  electric  field  can  lead  to  a 
considerable  increase  or  decrease  of  equilibrium  proton  charge  carrier  concentraticffis,  depending 
upon  electric  field  strength,  boundary  conditions  and  specimen  thickness. 

We  will  start  with  a  description  of  so-called  recombination  injection.  This  phenomenon  was 
originally  discovered  in  pure  water  (Petrenko  and  Qiesnakov  1990d)  and  later  s’lso  in  pure  ice 
(PetrerUco  and  Chesnakov  1990c).  A  reversible  increase  in  ac  conductivity  of  ice  and  water  thin 
layers,  through  which  a  weeik  direct  current  was  passed,  was  observed.  An  example  of  this  ef¬ 
fect  of  direct  current  on  water  conductivity  is  shown  in  Figure  26. 


Figure  26.  Dependence  ofac  conductivity  of  pure  water  a  on  time 
t,  when  a  small  dc  bias  (3  V)  ums  applied;  ac  conductivity  was  mea¬ 
sured  at  f=l  kHz;  T  =  20°C;  space  between  two  Pt  electrodes  was  3  mm. 
Dashed  lines  indicate  moments  when  the  bias  was  turned  on(t  =  0s)  and 
turned  off(t  s210s)  (after  Petrenko  and  Chesnakov  1990d). 


47 


Figure  27.  Alternating  current  conduc¬ 
tivity  oftuater(f=l  KHz)  nomudized  on 
equiliMum  conductivity  Oq  in  the  space 
between  two  electrodes.  TfK  upper  curve  is 
for  a  couple  afPt  electrodes,  while  the  lower 
curve  is  far  a  Pd  cathode  and  a  Pt  anode;  5- 
V dc  bias  is  applied;  T  =  20‘C  (after  Petrenko 
and  Chesrwkoo  1990d). 


The  condiictivity  increase  is  completely  reversible,  Le.,  it  ceases  after  the  dc  bias  is  switched 
off  and  is  not  related  to  electrodes  dissolving  in  water  or  ice.  The  eHect  is  equally  successfully 
observed  when  electrodes  made  from  Pt,  Au,  C  or  stainless  steel  are  used.  The  threshold  volt¬ 
age,  at  which  an  increase  of  a  arises,  appeared  to  be  equal  to  the  threshold  electrolysis  voltage 
(about  1.5  V  in  the  case  of  water  and  sli^tly  less  frian  2  V  in  ice).  It  has  been  shown  by  measur¬ 
ing  conductivity  between  two  micro-electrodes  that  the  conductivity  increase  is  not  uniform 
throughout  the  space  between  electrodes  (Fig.  27).  The  maximal  increase  of  a  is  observed  close 
to  the  anode  and  cathode,  damping  exponentially  as  distance  from  the  electrodes  increases.  To 
this  end,  the  characteristic  width  of  excess  conductivity  zones  coincides  well  with  the  difrusion 
length  of  ambipolar  H30'^  and  OH"  ions'  diffusion 


X.d=V2^  (162) 

where  D  is  the  diffusion  coefficient  of  the  less  mobile  ion  (OH")  and  Xi  is  the  non-equilibrium 
ions'  lifetime.  It  has  been  found  that  in  water  X.o=  3.4  x  10"*  m  (T = 20®C),  and  in  ice  4.5  x 
10"®  m  (T = -30®C).  Under  the  same  conditions,  the  non-equilibrium  ions'  lifetime  is  Tl = 7  sec¬ 
onds  in  the  case  of  water  and  Tl = 2  seconds  in  the  case  of  ice.  It  turned  out  that  an  increase  in  o 
in  zones  adjacent  to  the  electrodes  is  caused  by  formation  of  excessive  OH"  and  H30^  ions, 
concentrations  of  which  must  be  approximately  equal  owing  to  constraints  imposed  by  space 
charges. 

The  key  to  understanding  the  mechanism  of  excess  conductivity  was  provided  by  experi¬ 
ments  in  which  a  Pd  cathode  was  used.  Pd  absorbs  atomic  hydrogen,  generated  in  the  electrolysis 
process  according  to  the  reaction  defined  by  eq  23,  from  water.  In  cases  when  the  atomic  hy¬ 
drogen  does  not  get  into  water  (or  ice),  the  zone  of  excessive  conductivity  does  not  form  near 
the  cathode  (Fig.  27).  It  is  assumed  that  the  excess  ions  form  because  of  energy  release  in  the 
atomic  hydrogen  and  oxygen  recombination  process 


H -H  H -»  H2 -t- 4.5  eV 

(163) 

0-H0-»02-t-5.1eV. 

(164) 

The  quantities  43  and  5.1  eV  cire  adopted  from  the  book  of  Sokolov  (1976). 
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Figure  28.  Variations  in  l<m-frequenqf 
conductivity  of  pure  ice  measured  at  30 
Hzwhen  10-Vdc  bias  luas  turned  on  artd 
off;  7  =  -JO'C;  sample  thickness  is  0.5  mm; 
ohmic  electrodes  were  used  (after  Petrenko 
and  Chesnakov  1990c). 


Hie  energy  released  in  each  recombinaticHi,  defined  by  eq  163  and  164,  exceeds  by  many  times 
the  energy  required  for  producing  an  ion  pair  (~1  eV  in  the  case  of  ice  [see  the  Review  of  Experi¬ 
mental  Results  on  Ice  Conductivity  and  Dielectric  Permittivity  section]  and  059  eV  in  the  case  of  water 
[Hetcherl970]). 

Therefore,  we  can  assume  that  at  least  part  of  this  "tremendous"  energy  is  spent  on  produc¬ 
ing  new  ion  pairs.  Non-equilibrium  ions,  originating  near  the  anode  and  cathode,  then  diffuse 
into  the  specimen  bulk,  thus  increasing  the  bulk  conductivity.  For  this  phenomenon  to  be  ob¬ 
served  in  ice,  ohmic  electrodes,  described  in  the  Charge  Exchange  at  Ice/Metal  Interfaces  section, 
must  be  used.  The  effect  of  weak  electric  current  (-10^  A/cm^)  on  ice  conductivity  is  represented 
in  Figure  28.  Conductivity  was  measured  in  this  experiment  at  a  ftequency  of  30  Hz  « (%;  i.e., 
the  measured  conductivity  corresponds  to  static  conductivity.  The  latter  is  determined  by  OH“ 
and  H3<y'  ion  motion.  The  relative  change  of  high-frequency  conductivity  o„  is  practically  unde¬ 
tectable  in  this  case,  since  ~  Aa^. 

The  phenomenon  described  here  is  of  great  interest,  because  it  allows  us  to  determine  the 
lifetimes  and  diffusional  lengths  of  charge  carriers  in  water  and  ice  by  direct  methods. 

It  is  curious  that  for  thin  ice  (and  water)  specimens  having  thickness  L  less  than  diffusional 
length  A.D,  nonlinear  current-voltage  characteristics  are  observed  (Fig.  29).  They  can  be  readily 
explained  theoretically.  Let  us  consider  an  ice  specimen  with  L  «  such  that  ionic  concen¬ 
tration  «i  can  be  regarded  as  uniform,  independent  of  coordinates.  Also  let  a  be  the  number  of 
ion  pairs  produced  as  a  consequence  of  passing  a  single  electron  charge  e  through  the  specimen 
(a  >  1).  Then  the  balance  between  production  and  recombination  of  ions  will  be  defined  by  the 
equation 

L  (g  -pnj^)  -H  a  1  =  0  (165) 

where  G  is  the  rate  of  thermal  production  of  ions  in  ice  and  P  n  f-  (P  =  const)  is  the  recombina¬ 
tion  rate  in  the  bulk.  For  current  density  /,  the  following  relationship  holds 


'  =  (Iy)  ("l  +<>2)  J  =  {^)  [M  +1-21  "i  ^  ■ 

F  Ti  eq  165  and  166  we  readily  obtain  (for  weak  currents  a  ^  <  <  G ) 

i.e..  Ohm's  law.  For  large  currents  |a  ^  >  >  G  |  we  have 

/  =  a  oc 


(166) 


(167) 


(168) 
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Figure  29  illustrates  perfectly  the  trarisition 
frran  CSim's  law,  in  die  case  of  small  currents, 
to  square  current  dependence  from  voltage  as 
current  increases. 

In  concluscn,  let  us  note  that  the  energy  re¬ 
lease  near  the  cadrode  and  anode,  which  ob¬ 
viously  results  in  an  increase  in  die  number  of 
ions,  may  also  increase  die  number  of  ^errum 
defects  in  comparable  quantities.  However, 
since  the  concentration  of  the  latter  exceeds  by 
several  orders  of  magnitude  die  ions'  ccxioen- 
traticxi  in  pure  ice,  diis  does  not  alter  hb  appre¬ 
ciably.  It  appears  that  for  Bjerrum  defects  the 
phenomenon  described  next  is  mudi  more  es¬ 
sential. 


RECOMBINATION  EXTRACTION  OF 
CHARGE  CARRIERS  FROM  ICE 


Figure  29.  Current-voltage  characteristic  of  a 
thin  (20-fm)  ice  specimen  with  ohmic  electrodes; 
T  =  -30‘C;  electrode  area  2  cn^  (after  Petrenko 
and  Chesnakcv  1990c). 
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First,  a  decrease  in  a  static  electric  field 

being  applied,  was  observed  by  Petrenko  et  al. 

(1983)  (see  Fig.  2  of  their  paper).  However,  at 

that  time  no  importance  was  attached  to  this  observation.  Later,  an  efiect  of  static  electric  fields 
CHI  ice  ccHiductivity,  which  is  opposite  to  that  just  described  in  the  previous  section,  was  observed 
and  studied  in  th<»  work  of  Petrenko  and  Schulson  (1992a,  b).  It  appeared  that  high-frequency 
conductivity  (co  >  Wd)  decreased  by  a  factor  between  2  and  10  while  a  relatively  small  direct  volt¬ 
age  was  being  applied  to  pure  ice  specimens  (Fig.  30).  Hie  relative  decrease  of  when  a  di¬ 
rect  electric  field  is  applied  becomes  noticeable  starting  at  a  temperature  of  -10°C,  and  it  inctease 
rapidly  at  lower  temperatures. 

The  phenomenon  is  practically 
insensitive  to  the  type  of  elec¬ 
trodes  used.  So,  it  has,  to  the  first 
approximation,  the  same  magni¬ 
tude  both  for  ohmic  electrodes  as 
well  as  for  common  stainless 
steel  electrcxles;  i.e.,  the  effective 
charge  exchange  between  ice  and 
electrodes  does  not  play  such  a 
significant  role  as  it  did  in  the 
case  of  recombination  injection. 

This  observation  is  well  illus¬ 
trated  by  Figure  30.  As  can  be 
seen  from  this  figure,  the  value  of 
the  observed  decrease  in  is 
comparable  for  both  dc  bias  po- 

larities.  At  the  same  time,  dc  ftSxtcdO.  Depaidsm^si^mmMlmdhi^feq^comluc- 

flows  leadUy  it.  the  positive  di-  “"f  f  “1“^  '  T  T  f  A'  , 

^  Single  crystal  of  ice:^  =  -40°C;  stainless  Steel  ohmic  electrodes  were 

rection,  when  the  electrodes  op-  p,.  positive  and  negative  polarities  cf  current  vary  in  20 

erate  in  the  ohmic  mode,  and  (after  Petrenko  and  Schulson  1992a,  b). 
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practically  does  not  flow  in  tfw  c^^xjsite  di¬ 
rection.  In  terms  of  potential  distributum, 
this  means  fliat  in  flie  first  case  the  potential 
difference  is  applied  to  the  ice  bulk  and  in 
fire  second  case  it  almost  completely  drops 
at  two  ice/metal  interflaxs. 

Ihe  dependence  of  ice  conductivity  mea¬ 
sured  at  high  frequency  on  specimen  thick¬ 
ness  is  shown  in  Hgure  31.  Ohmic  electrodes 
have  been  employed  in  fitis  experiment  and 
an  electric  field  has  been  applied  in  a 
"straig^' directkm.  An  irKiease  in  a;,  for  ice 
thickness  less  than  10-20  pmiscausedby 
the  efrect  of  recombination  ir^ection  (which 
predicts  o„  ~  Lr^,  see  the  previous  section, 
eq  168),  and  also  by  the  effect  of  rou^uiess 
present  at  the  surface  of  the  electrodes.  As 
for  file  characteristic  a..  irKrease  in  the  speci¬ 
men  thickness  range  from  25  to  100  pm,  it  is 
easily  reproducible  and  evidently  reflects 
the  internal  distribution  of  majority  charge  carriers  in  a  thin  ice  specimen.  To  account  for  the 
dependerve  oJ^L)  shown  in  Figure  31,  the  authors  of  this  study  assumed  that  ke/metal  inter¬ 
faces  are  "drains"  for  Bjerrum  defects;  i.e.,  the  recombinati(H\  rate  of  L~  and  defects  at  the 
interfaces  is  much  hig^  than  in  the  bulk.  As  a  result,  the  charge  carrier  concentrations  decrease 
profoundly  in  layers  in  which  thickness  is  comparable  to  diffusion  lengths.  When  the  potential 
difference  is  applied,  the  recombination  rate  at  flie  interfaces  increases  even  more  for  two  rea¬ 
sons.  The  first  is  the  increase  of  defect  densities  npand  wl  inside  layers  of  space  charge  of  thick¬ 
ness  about  the  screening  length  S 1  pm,  since  the  rate  of  recombinaticxi  is  proportional  to  («d”l)- 
The  second  is  additicmal  drift  of  carriers  in  the  electric  field  to  file  interfoces.  Simple  numerical 
evaluaticHis  performed  by  Petrenko  and  Schulscm  may  confirm  the  assumed  mechanism.  Al¬ 
though  foe  exact  mechanism  of  flie  direct  electric  field  effect  on  is  not  absolutely  dear  yet, 
the  fact  of  essential  and  reversible  decrease  in  majority  charge  carrier  concentrations  is  reliably 
established.  This  phenomenon  promises  interesting  applicatiOTis  in  studies  of  ice  properties  that 
depend  upon  majority  carrier  concentrations — foe  plastic  flow  and  inelastic  relaxation  of  ice, 
for  instance. 


Figure 31.  Dependence c^a^on  sample  thickness  L; 
T = -40X1;  single  crystalline,  pure  ice  samples  used 
with  dimic  electrodes  (after  Petrenko  and  Sdtulson 
1992a). 


PROTON  INJECTION  FROM  Pd  ELECTRODES  INTO  ICE 

The  method  of  injection  current  measurement  has  been  very  effidently  employed  for  a  long 
time  in  foe  study  of  electrical  properties  of  ordinary  semicondudors  and  dielectrics  (Lampert 
and  Mark  1970).  By  studying  foe  charge  carriers'  ir^ection  into  foe  bulk  of  the  investigated 
material,  one  can  determine  such  important  parameters  as  charge  carrier  mobility,  foeir  equi¬ 
librium  concentrations,  foe  energy  position  of  traps  and  foe  traps'  concentrations. 

In  this  method  the  carriers'  space  charge  is  "injected"  into  a  material  from  an  electrode.  The 
charge  value,  and  hence  foe  number  of  charge  carriers,  is  known  because  those  quantities  are 
determined  by  the  laws  of  electrostatics.  Since  foe  concentraticm  is  already  known,  then  foe 
mobility  can  be  calculated  from  conductivity  measurements. 

Let  us  clarify  foe  main  ideas  of  such  experiments  by  following  simple  reasoning.  Suppose  a 
direct  electric  voltage  V  is  applied  to  a  conductor  with  two  electrodes  (Fig.  32).  The  electrode 
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Figure  32.  Injection  of  protons  into  ice  from 
Pd  electrode  saturated  by  hydrogen.  Pd  cath¬ 
ode  works  as  proton  acceptor.  Upper  diagram  il¬ 
lustrates  distribution  cf  protonic  concentration 
n  and  electric  field  E  inside  the  specimen  bulk. 


under  positive  potential  is  able  to  inject  positive  charge  carriers  with  charge  value  e;  i.e.,  this 
electrode  is  an  "infinite  capacity"  reservoir  of  sudr  carriers.  Therefore,  the  carrier  concentration 
in  the  vicinity  of  the  electrode  m(+0)  =  «>,  which  implies  that  the  electric  field  E(+0)  =  0,  because 
of  "infinite"  conductivity  of  the  media  at  this  point.  The  other  electrode,  at  the  point  x  =  L,  is 
the  charge  carrier  receiver.  Apparently,  when  a  voltage  V  is  applied,  foe  carriers  will  be  "pushed" 
into  foe  bulk,  increasing  the  conductivity  of  a  specimen.  In  foe  case  when  concentrations  of  in¬ 
jected  carriers  are  much  greater  than  equilibrium  ccmcentrations,  we  can  readily  write  down  a 
closed  system  of  equations  that  defines  injection  current  density  /,  electric  field  strength  £  and 
carrier  concentrations  n 


—  =  —  =  ^  (Poisson's  equation)  (169) 

dx  eeo  ^ 

f  =ej  =  enpE  =  const  (170) 


£dx  =  P;  £  (+0)  =  0;  «  (+0)  =  oo . 

'o 


(171) 


The  solution  is  given  by  an  elegant  formula 

]  =  (172) 

8  1,3 

in  which  there  appears  only  one  unknown  parameter,  p,  which  can  be  found  from  experimen¬ 
tal  measurement  of  /.  Therefore,  we  can  determine  charge  carrier  mobility  by  means  of  injec¬ 
tion  current  measurements.  That  is  not  all.  The  concentration  of  injected  carriers  b 

«  1  LJQ.v'Jl  (173) 

4  e  ^ 
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and  die  "average"  bulk  ccsicentratian  is 


L 

<n>  =1  f  =2.  i_!0  V  (174) 

0  ^  ^ 

If  equilibrium  diarge  earners  widi  concentration  rq  are  present  in  die  investigated  material,  dien 
at  small  voltages  V  (small  level  of  injecticHi)  a  linear  zone  will  be  observed  on  die  voltage-cur¬ 
rent  characteristic,  i.e.,  CMim's  law 


/o=enno^.  (175) 

la 

Hie  transition  from  dependence  (eq  175)  to  eq  172  will  occur  in  die  region  where  /  =  Jq.  From 
here  it  follows  that 


V^tran  =*  J  «0  (176) 

\9  e  tQ  I 

i.e.,  having  determined  from  current-voltage  characteristics,  we  can  calculate  Hq.  So,  the 
ii^ection  currents  technique  allows  us  to  determine  charge  carrier  concentrations  and  mobili¬ 
ties  in  diose  materials  in  which  the  Hall  effect  cannot  be  measured — because  of  extremely  small 
mobility,  for  instance,  as  in  the  case  of  ice. 

At  first  glance  die  injection  of  protonic  charge  carriers  into  ice  is  extremely  complicated,  since 
there  are  four  types  of  such  carriers  in  ice  and  their  motion  is  correlated  through  a  configura¬ 
tion  vector  In  fact,  the  madiematical  description  of  injection  currents  theory  in  ice  is  compli¬ 
cated  and  space-consuming.  The  results  of  this  dieory  are  nevertheless  simple  and  dear,  which 


allows  us  to  interpret  them  successfully  in  the 
context  of  the  consideration  just  performed. 

The  theory  of  ir^ection  currents  in  ice  was  de¬ 
veloped  for  the  first  time  in  the  work  of  Petrenko 
et  al.  (1983),  in  which  a  solution  was  obtained  for 
one  particular  case:  the  injected  proton  concentra¬ 
tion  exceeded  the  equilibrium  ion  concentration, 
but  was  still  much  less  than  die  Bjerrum  defect 
concentration.  The  general  problem  was  solved 
later  (Petrenko  and  Ryzhkin  1984b).  Either  a  pro¬ 
ton  can  be  ir^ected  into  ice  fi:om  outside,  which 
introduces  an  H30‘''  ion  and  a  D-defect  si¬ 
multaneously  (see  Fig.  8),  or  protonic  "hole"  oc¬ 
curs,  which  is  equivalent  to  the  production  of  an¬ 
other  defect  pair — an  OH"  ion  and  an  L-defect. 
The  equations  describing  the  protexi  and  protonic 
"hole"  irqection  into  ice  are  the  same.  Therefore, 
it  is  enough  to  consider  the  case  of  proton  injec¬ 
tion  as  an  example. 

The  theoretical  current-voltage  characteristic 
of  injection  protonic  currents  in  ice,  adopted  from 
the  paper  of  Petrenko  and  Ryzhkin  (1984b),  is 
shown  in  Figure  33.  Zone  I  (Ohm's  law)  is  ob¬ 
served  when  die  average  concentration  of  ir^ected 


Figure  33.  Theoretical  current-voltage  charac¬ 
teristic for  proton  injection  into  ice  (see  expla¬ 
nations  in  the  text  [after  Petrenko  and  Ryzhkin 
1984b]). 
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protons  Ftp  is  much  less  than  equilibrium  kmk  coixsentFatkin  Hj.  VVhen  Kp  becomes  gireato' than 
Fii,  but  is  ^  less  rium  equilibrium  Bjenxim  defect  ctmcentration  Fig,  zone  n  must  be  obsoived, 
in  which  /  ~ 


J  -  Egeom 
Is  e^L^I 


(177) 


As  can  be  seen  horn  tttis  fonnula,  we  can  detennine  H3Q*’  ion  mobility  from  measurement 
of  /.  Note  also  riiat  the  magnitude  of  injection  currents  in  zcme  n  is  proportional  to  a  large 
(- 10^)  static  dielectric  permittivity  of  ice  e^.  The  magnitude  of  the  transition  voltage  Vfi  is  char¬ 
acterized  by  ttie  condition  diat  Fip  »  ff^.  Hence,  using  eq  176  and  e  =  e,,  one  can  determine  the 
equilibrium  ion  concentraticHi. 

When  voltage  magnitude  and  level  of  injectkm  are  increased,  it  may  happen  that  rip  ^  Fig-  In 
this  case,  since  for  ii^ected  carriers  fii  =  rig  holds,  dte  concept  of  majority  and  minority  charge 
carriers  ceases  to  exist.  The  constraint  (m  carrier  moticai  imposed  by  the  configuration  vector  is 
also  eliminated,  die  carrier  motion  becomes  uncorrelated,  and  e  drops  from  e,  » 100  to  e,.  a  32. 
This  leads  to  the  following  observation  in  zone  IV,  where  the  current  magnitude  is  die  square 
of  the  voltage,  but  the  proportionality  coefficient  is  smaller 


/  =• 


yeini+eajia)  ^3 


(178) 


This  allows  us  to  calculate  both  the  minimal  mobility  of  the  two  (pi  or  P3)  fiom  measurement 
of  /  and  the  equilibrium  ^errum  defect  concentraticn  from  the  value. 

If  measurements  of  current  in  die  case  of  protonic  '"hole"  injection  were  performed,  we  could 
determine  all  four  mobilities  Pi,...,m  and  aU  four  equUibrium  concentrations  However, 

experimental  studies  in  the  field  of  ir^ection  currents  are  much  less  successful  than  theoretical 
aies.  In  all  experimental  works  on  ir^ection  currents  in  ice,  hydrogen-saturated  Pd  or  palladium 
black  were  u^  as  proton  injecting  electrodes.  Here,  hydrogen,  which  could  be  dissolved  in 
Pd  with  concentrations  close  to  one,  gives  its  single  electron  to  Pd  and  readily  migrates  in  the 
Pd  lattice.  Pd  can  be  saturated  with  hydrogen  in  a  water  electrolysis  process  (when  Pd  is  a  cath¬ 
ode)  or  in  a  heating  process  imder  high  pressure  in  a  hydrogen  atmosphere.  Such  protonic  in- 


Figure  34.  Dispersion  curves  for 
pure,  monocrysialline  ice  speci¬ 
men  with  two  H-loaded  Pd  elec¬ 
trodes  of  area  =  15  x  9  mm^, 
thickness  =  0.7  mm,  tempera¬ 
ture  =  -26‘C.  C(0)  and  G(0)  are 
the  specimen  capacitance  and  con¬ 
ductance  without  dc  bias;  C(V) 
and  G(y)  are  the  same  with  dc  bias 
ofV  =  2.0  kV  (after  Petrenko  et  al. 
1983). 
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jecticot  into  ice  was  first  carried  out  by  Engelhart  and  Riehl  (1%5, 1966),  who  measured  current- 
voltage  characteristics  of  thin  (-  05  ii\m)  ice  specimens  with  Pd  electrodes.  In  this  study  they 
were  not  able  to  obtain  zrmes  in  which  current  was  distinctly  proportional  to  the  square  of  the 
voltage,  and  the  current-voltage  characteristics  were  obviously  complicated  by  the  effect  of 
proton  traps  present  in  the  imperfect  ice  crystals  used  in  those  experiments.  Neverti\eless, 
quantitative  evaluations  that  can  be  made  on  the  basis  of  existing  theories  indicate  fiiat  in  those 
first  studies,  the  irqection  of  protcMis  into  the  bulk  was  undoubtedly  achieved.  The  same  is  not 
true  for  titeir  later  research  (Engelhart  et  al.  1%9),  performed  on  ice  specimens  of  such  large 
thickness  (10  cm)  that,  although  the  zone  of  square  dependence  was  observed  in  current-volt¬ 
age  characteristics,  it  cannot  be  explained  in  terms  of  space  proton  injection  (/  ~  Lr^). 

To  prove  that  ice  bulk  conductivity  (Oj  and  o*)  increases  as  a  result  of  proton  injection, 
Petrenko  et  aL  (1983)  made  measuremaits  at  ac  along  with  dc  measurements  tluit  elirrurrate  any 
doubts  concerning  the  effect  of  ice/metal  interface  resistance  on  current-voltage  characteristics. 
The  changes  in  a({0)  and  /  due  to  proton  injection  into  ice  are  shown  in  Figures  34  and  35.  Us¬ 
ing  stationary  currents  of  monopolar  proton  injection  made  it  possible  for  H30‘*’  ion  mobility 
to  be  determined  (see  Table  2). 

A  remarkable  experimental  study  on  the  measurement  of  transient  injection  currents  in  ice 
was  carried  out  by  Eckener  et  al.  (1973).  In  these  experiments  a  short  laser  pulse  released  pro¬ 
tons  from  a  thin  palladium  black  film,  which  was  kept  in  ccmtact  with  ice.  The  released  protons, 
under  the  action  of  force  from  a  constant  electric  field,  were  moving  from  anode  to  cathode,  thus 
producing  a  current  pulse,  the  duration  of  which  was  equal  to  the  transit  time  of  protons  through 
an  ice  specimen 

t  =  _L_  =-LL.  (179) 

niE  mV 

The  value  of  H30^  ion  mobility,  which  is  obviously  file  most  reliable  data  for  today,  has  been 
obtained  by  this  method  (see  Table  2). 

The  extensive  use  of  the  injection  current  technique  in  ice  is  limited  by  basic  difficulties  in 
preparation  of  electrodes  that  would  efficiently  inject  protons  into  ice.  Contamination  processes 
at  the  ice/Pd  interface,  which  are  poorly  controlled,  deteriorate  the  reproducibility  of  results 
considerably. 


Figure  35.  Direct  current  I,  ca¬ 
pacitance  C  and  conductance  G 
at  100  Hz,  as  functions  of  bias 
voltage  V  for  specimen  with  two 
H-loaded  Pd  electrodes  of  area  = 
6  x9  mm^,  thickness  =  0.4  mm, 
temperature  =  -20°C.  Broken 
lines  of  slopes  1  and  2  have  been 
drawn  for  the  current  characteris¬ 
tics  (after  Petrenko  et  al.  1983). 
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HELD  EFFECT  TRANSISTOR  MADE  OF  ICE 


The  operation  principle  of  the  field  effect  transistor,  if  one  tries  to  describe  it  in  a  few  words, 
is  as  follows.  Let  us  consider  a  conductor/dielectric  interface,  as  shown  in  Figure  36.  Suppose 
that  normal  to  this  interface,  a  static  electric  field  is  applied,  £(-0)  =  which  attracts  charge 
carriers  to  the  interface.  A  thin  layer  of  excess  concentration  Nj  (dimension  m“^)  will  form  in 
the  conductor  along  the  interface.  Since  the  electric  field  at  the  infinity  £(+<»)  =  0,  then  fi'om  the 
well-known  Gauss  electrostatic  theorem  it  follows  that 

£1  =  (180) 

EdCO 

where  is  the  dielectric  pjermittivity  of  the  dielectric.  Then  the  surface  conductivity  of  the  con¬ 
ductor,  resulting  from  the  applied  field,  is 

Xs  =eiNmi  =£deo£xFi-  (181) 


Formula  181  shows  how  the  value  of  charge  carrier  mobility  Hi  can  be  determined,  if  the  ap¬ 
plied  field  £j^  is  known  from  the  experimental  conditions  and  A,  is  measiued.  If  there  are  two 
types  of  charge  carrier  present  in  the  conductor  (for  instance,  H3O*  and  OH“  ions  in  ice),  fiien 
in  order  for  the  second  type  of  charge  carrier  mobility  to  be  determined,  we  should  just  change 
the  polarity  of  EJ 

In  case  of  ice,  this  apparent  simplicity  is  complicated  by  the  presence  of  two  more  charge 
carriers;  D-  and  L-defects.  However,  the  effect  of  the  electric  field  on  their  concentration  is  neg¬ 
ligibly  small  compared  to  the  change  in  ion  concentrations.  This  fact  is  easy  to  understand  if 
we  recall  that  the  carrier  concentration  ni(x)  depends  upon  ffie  magnitude  of  electric  potential 
at  a  given  point  (x),  as 


n; 


«, exp 


(  ei(p[x)\ 

1  ksT  I 


(182) 


where  <p(x)  is  the  electric  potential  and  is  the  carrier  concentration  at  infinity,  where  there  is 
no  electric  field  and  cp  =  0. 


Since  the  value  of  the  ions' 
electric  charge  Ci  =  0.62  e  is 
much  greater  than  0.38  e, 
which  is  the  Bjerrum  defect 
charge,  then  the  effect  of  the 
electric  field  is  also  much 
more  profoimd  for  ions. 

As  long  as  Ojon  <  Oq  in  the 
surface  layer,  the  static  con¬ 
ductivity  will  be  determined 
by  ions,  so  that 

Xs  =  (^)2e^eotii£x  (183) 
where  the  factor  (e/ei'^  ~  2.6 


follows  firom  eq  70,  which  d  -  Figure  36.  Dieleciric/conductor  interface  with  applied  static  electric 
scribes  the  static  conductivity  fjgj(j  77, ^  diagram  below  shows  schematically  a  distribution  of  charge 

of  ice.  When  E^_  is  increasing,  carrier  concentration  n. 
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Oio„  can  exceed  since  the  ionic  con¬ 
centration  grows  faster  than  Bjerrum 
def^  concentration.  Bjerrum  defects 
cease  to  be  majority  charge  carriers  and 
their  mobility  and  charge  appear  in  eq 
183.  So,  by  means  of  measuring  Og  at 
different  magnitudes  and  polarities  of 
electric  field  E^,  we  can  determine  all 
four  mobilities. 

The  field  effect  transistor  made  of 
pure  ice  was  experimentally  realized  by 
Petrenko  and  Maeno  (1987).  The  layout 
of  such  a  transistor  is  shown  in  Figure 
37.  A  thin  (about  1000  A)  dielectric  layer 
of  SiQz  is  placed  onto  a  conducting  sub¬ 
strate  made  of  Si.  Gold  r’ectrodes  are 
deposited  onto  the  Si02  layer  by  means 
of  photolithography,  and  are  then  used 
for  surface  conductivity  measurements. 

To  increase  tne  total  length  of  contact 
and  thus  reduce  the  measured  surface 
conductivity,  the  electrodes  had  the 
shape  of  two  combs  inserted  tooth- 
wards  into  each  other,  as  shown  in  the 

upper  section  of  Figure  37.  A  bias  was  applied  between  the  Si  substrate  and  ice,  and  a  layer  was 
formed  at  the  SiOj/ice  interface  with  an  excess  concentration  of  either  H3CP'  ions  ("plus"  is  at 
ice)  or  OH"  ions  ("plus"  is  at  Si). 

Figure  38  illustrates  the  alternation  of  the  surface  conductivity  of  ice,  when  rectangular  volt¬ 
age  pulses  of  different  polarity  are  applied.  The  values  of  and  P2  determined  in  this  experi¬ 
ment  can  be  found  in  Table  2. 


mm 

Figure  37.  Field  ^ct  ice  transistor  (after  Petrenko  and 
Maeno  1987). 
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Figure  38.  Time  diagrams  of  ice  surface  conduc¬ 
tivity  and  dc  bias  (V)  applied  to  the  transistor 
shown  in  Figure  37;  T  =  -33.1  “C;  surface  conduc¬ 
tivity  was  measured  at(=WHz  (after  Petrenko  and 
Maeno  1987). 
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"CROSSOVER"  IN  THE  DIELECTRIC  PERMITTIVITY  OF  ICE 

Until  now,  we  have  referred  to  static  dielectric  permittivity  of  ice  as  a  large  quantity,  about 
10^.  For  the  most  part  this  is  true.  But  in  some  instances,  at  particular  temperatures  and  impu¬ 
rity  concentrations,  can  decrease  drastically  and  achieve  small  values,  comparable  with  those 
of  =  3.2.  This  phenomenon,  which  has  become  known  as  "crossover,"  can  be  successfully 
explained  in  the  context  of  Jaccard's  model  of  the  electrical  conduchvity  of  ice. 

If  we  examine  attentively  formula  72,  which  describes  as  a  function  of  temperature  ((t>  = 
3.85  k]gT roo)  and  conductivities  Oj^n,  then  we  will  nohce  that  if  the  following  relationship  is  sat¬ 
isfied 


Oi  4-  02  _  03  +04  (184) 

ei  C3  ^ 

the  second  term  in  formula  72  reduces  to  0  making  equal  to 

The  physical  interpretation  of  such  a  significant  drop  in  Cg  is  that  if  the  above  condition  (eq 
184)  is  satisfied,  then  the  charge  carriers  moving  under  the  action  of  an  external  electric  field — 
the  fluxes  of  ions  and  Bjerrum  defects — ^appear  to  be  equal  (see  eq  64).  Therefore,  the  configu¬ 
ration  vector  Q  =  0  all  the  time,  i.e.,  the  predominant  orientation  of  hydrogen  bonds  is  not  formed. 
Consequently,  there  is  no  additional  returning  force  from  D  acting  upon  the  charge  carriers.  It 
is  this  force,  resulting  from  D,  that  converts  pairs  of  major  charge  carriers  into  "elastic  dipoles." 
In  other  words,  under  the  condition  of  eq  183,  inasmuch  as  Q  =  0,  the  correlation  in  charge  car¬ 
rier  motion  ceases  and  they  will  move  in  the  same  fashion  as  in  an  ordinary  conductor.  In  this 
case,  no  further  contribution  to  apart  from  (dielectric  permittivity  of  media)  can  be  expected! 

In  very  pure  ice  the  "crossover"  phenomenon  is  not  observed,  since 

Oion  =  (Oi  +  (^)  <  <  (03  +  04)  =  OB  .  (185) 

This  is  related  to  the  fact  that  at  high  temperatures  (0  to  -10®C),  Ob  =  exceeds  by  about  two 
or  three  orders  of  magnitude  the  ionic  conductivity  Ojon  =  Og.  If  the  temperature  is  de¬ 
creased,  Og  and  Ojon/  since  they  have  comparable  values  of  activation  energies — about  0.5-0.6 
eV — will  never  become  equal  until  they  reach  the  lower  limit  of  ac  measurements  possible  in 
pure  ice  (=  150  K). 

However,  the  picture  changes  drastically  if  there  are  impurities  present  in  ice,  since  their  ef¬ 
fect  at  Cion  and  Ob  is  different  (see  the  Concentration  of  Charge  Carriers  section). 

So,  for  example,  when  H30'^  "impurity"  ions  are  introduced,  is  almost  constant  at  high 

temperatures,  while  Ob  decreases  rapidly,  as  temperature  decreases.  Temperature  dependen¬ 
cies  for  Ob  and  Ojon  in  ice  doped  with  a  small  amoimt  of  HF  (1 .3  x  10^^  cm"^)  are  shown  in  Fig¬ 
ure  39.  Close  to  the  point  T„  =  230  K,  Ob  and  tjjon  become  equal,  and  at  lower  temperatures  they 
change  places:  Ob  <  0,00.  In  other  words,  at  the  point  of  intersection  ("crossover"),  the  ions  turn 
from  minority  into  majority  charge  carriers,  and  Bjerrum  defects  become  the  minority  charge 
carriers.  So,  we  have 


at  T  >  T„:  o„=  Ob  ;  <%  = 

(186) 

at  T  <  ,  Ooo  “  qion  ;  Og  =  OB- 

(187) 

Figure  40  illustrates  a  sharp  decrease  of  eg  in  the  vicinity  of  Tq.  in  HF-doped  ice.  Hubmann  (1978) 
observed  "crossover"  in  NHj-doped  ice  (5  x  10"®  mol/L)  when  the  pressure  was  increased,  since 
the  latter  increased  Ojon  and  decreased  Ob- 
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Figure  39.  Temperature  depen¬ 
dences  ofH^  ion  conductivity 
(T,  and  Bjerrum  defects'  conduc- 
tivity  Opi  calculated  from  experi¬ 
mental  results,  using  best-fit 
method.  Ice  doped  by  HF  in  the 
concentration  of  (1.3  ±  0.1)  x  id^ 
molecules/n^  (after  Complin  et  al. 
1978). 


m  “Pure' Ice 

•  (1.1+ 0.2  )x10-«  Mole/Ll 
^  (8.0^0.2)x10-6Mole/L: 
^  (5.1  ^0.2)  X  10-5  Mole/L- 
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Figure  40.  A£  =  Eg  -  £^as  a  func¬ 
tion  of  temperature  for  ice  doped 
with  HF  concentrations  (after 
Cantplin  and  Glen  1973.) 


A  remarkable  example  of  "double  crossover"  can  be  found  in  an  old  report  by  Steinemann 
(1957) — ^which  became  a  classic — ^who  studied  dielectric  permittivity  of  HF-doped  ice.  At  the 
high  temperatures  used  by  Steinemann  (-3®C),  it  is  possible  to  pass  through  the  "crossover" 
point  twice,  since  Og  changes  by  an  amount  proportional  to  wnp/  and  Gjon  by  an  amount  pro¬ 
portional  to  VriHF-  Therefore,  the  equation 

OB  =  <^on-»  OBO+'*1«HF  =  <^o  + 
can  be  satisfied  at  two  different  values  of  (see  Hg.  41). 


Concentration  n  hf  of  HF  (mol  I  “') 


Figure  41.  Static  relative  permittivity  and  dielectric  relaxation  time  t 
at-WXIasa  function  of  the  concentration  n^F  in  the  solution  from  which 
the  polycrystalline  ice  was  formed  (after  Steinemann  1957). 

THERMALLY  STIMULATED  DEPOLARIZATION 

The  method  of  Thennally  Stimulated  Depolarization  (TSD)  is  successfully  used  in  studies 
of  dielectrics  and  materials  with  low  conductivity.  Its  generally  accepted  advantages  are  top 
instrumental  simplicity  and  high  sensitivity.  Figiue  42  illustrates  the  basic  principles  of  TSD 
measurements. 

In  this  method  the  material  under  investigation  is  first  cooled  imder  conditions  of  an^exter- 
nal  electric  field  E  being  applied.  Under  the  acficari  of  this  field,  an  electrical  polarization  P  arises 
in  the  material  (ice  in  our  case),  which,  after  cooling  to  a  sufficiently  low  temperature,  appears 
to  be  "frozen,"  i.e.,  it  remains  after  the  external  field  is  switched  off.  Here,  the  concept  of  "low 
temperature"  is  relative,  since  for  alkaline  halide  crystals  room  temperature  is  already  sufficiently 
low. 

When  the  specimen  is  heated  (usually  at  a  constant  rate  dT/dt  =  y)ti\e  frozen  electric  dipole 
moment  relaxes.  An  electric  current,  usually  measured  by  an  electrometer  and  by  an  XY  re¬ 
corder,  flows  through  the  external  circuit.  The  analysis  of  the  /(f)  recording  frequently  allows 
us  to  detemune  important  parameters  of  the  conductivity  processes. 

Let  us  illustrate  this  with  a  simple  example.  The  relaxation  rate  of  electric  polarization  P  can 
be  described  in  terms  of  simple  exponentuil  time  dependence 

P=Poexp{-tj.  (188) 

The  relaxation  time  t  as  a  rule  is  defined  in  terms  of  some  thennoactivation  process,  so  that 

T  =  To  exp  I  (189) 

IksT  I 

where  E^  is  an  activation  energy. 
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Figure  42.  Schemes  of  polarization 
(cooling)  and  depolarization  (heating) 
in  the  TSD  technique. 


As  we  have  already  seen  in_^  Relaxation  Times  of  the  ElecMc  Polarization  and  Electric  Fields  in 
Ice  section,  the  current  density  /  and  dielectric  polarization  P  are  related  through  the  simple  re¬ 
lationship  (eq  151) 

dt 


and,  taking  into  consideration  eq  188,  we  obtain 

Now,  substituting  eq  189  into  190  and  considering  that 


(190) 


P  =  -\j{t)dt 


(191) 


we  ultimately  obtain 


Jit) 

oo 

\j{t)<it 


(192) 


From  this  it  can  be  seen  that  if  a  graph  of  the  natural  logarithm  of  the  left  side  of  eq  192  is 
plotted  versus  {  i/IcqT  },  then  it  will  be  a  straight  line  witi\  slope  EfJ.  That  is,  in  this  way  the  value 
of  die  activation  energy  of  the  charge  transfer  process,  which  leads  to  relaxation  of  P,  can  be 
determined.  Furthermore,  since 


Po=\jit)clt 

o 

tiien  the  area  below  the  bell-shaped  curve  of  TSD  (Fig.  42)  contains  mformation  about  the  num- 
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ber  of  "polarized  and  frozen"  dipoles,  P  being  the  result  of  tfveir  reorientation.  Those  are  the  basic 
principles  underlying  the  TSD  method.  Quite  often  the  simple  TSD  analysis  presented  above 
is  complicated  in  practice  by  the  presence  of  several  relaxation  processes  instead  of  one.  For 
example,  the  initial  "conserved"  polarization  might  be  defined  in  terms  of  reoriortation  of 
several  different  corKentrations  and  activation  energies  of  elementary  dipoles. 

In  the  heating  process  different  dipoles  will  relax  at  difierent  temperatures.  Thus,  instead  of 
one  "simple"  bell-shaped  TSD  current  pulse,  we  will  observe  several.  For  the  analysis  of  such 
situations,  several  special  techniques  have  been  developed  that  allow  us  to  separate  one  relax¬ 
ation  process  from  another.  Let  us  look  at  some  of  those. 

If  it  is  possible  to  guess  by  the  shape  of  the  curve  /(f)  how  many  peaks  there  are  on  this  curve, 
then  we  can  employ  computer  analysis,  seeking  to  find  the  best  fitting  theoretical  curve. 

Another  method  that  allows  us  to  amplify  separate  "peaks"  on  the  background  of  other  peaks 
is  switching  on  the  external  field  £,  not  during  the  whole  cooling  period  of  the  investigated 
specimen,  but  only  in  that  temperature  range  in  whidt  ttiis  peak  will  "develop"  afterwards.  Then 
higher  temperature  peaks  will  not  appear,  because  fire  temperatures  are  too  low  for  their  pro¬ 
duction.  TSD  "peaks"  at  lower  temperatures  wiU  have  enough  time  to  relax  during  the  cooling 
process,  after  the  external  field  has  been  switched  off. 

A  modification  of  this  method  is  thawing  of  different  dipoles  systems  in  turn;  i.e.,  when  the 
temperature  is  increasing  and  after  the  first  (the  lowest  temperature)  peak  is  pcissed,  the  tem¬ 
perature  is  decreased  again  and  this  already  "developed"  peak  does  not  contribute  to  a  new 
thawing  cycle.  Then  this  procedure  is  repeated  until  the  next  peak,  and  so  on. 

For  ice  research  the  TSD  method  is  attractive  because  it  provides  a  way  to  study  conductiv¬ 
ity  processes  at  temperatures  so  low  that  ac  measurements  become  either  impossible  or  strongly 
impeded.  The  research  in  TSD  on  pure  and  doped-ice  monocrystals  started  over  two  and  a  half 
decades  ago  (see  EngeUiart  and  Riehl  1965,  Gelin  and  Stubbs  1965).  Detailed  experimental  studies 
of  TSD  in  ice  were  carried  out  by  Bishop  and  Glen  (1969),  Loria  et  al.  (1978),  Johari  and  Jones 
(1975)  and  Apekis  and  Pissis  (1987).  In  the  last  report  the  reader  can  fiind  an  extensive  list  of 
papers  concerned  with  TSD  in  ice.  In  many  papers  the  observation  of  several  TSD  peaks  has 
been  reported,  the  magnitude  and  location  of  which  depend  upon  doping  impurity  types  and 
concentrations  (see  Zaretskii  et  al.  1987b,  for  instance).  Two  examples  of  "thermograms," 
adopted  from  this  report,  are  shown  in  Figure  43.  In  the  reports  of  various  authors  concerned 
with  ice  TSD,  a  considerable  discrepancy  can  be  foimd  in  activation  parameter  values,  calcu¬ 
lated  from  those  "thermograms,"  and  their  various  interpretations.  We  deliberately  avoid  com¬ 
menting  on  these  numerical  data  for  the  reason  explained  below. 

In  all  the  dted  experimental  reports,  the  analysis  was  performed  on  the  basis  of  simple  model 
considerations  just  discussed.  However,  such  considerations  are  not  applicable  in  the  case  of 
ice,  because  of  the  presence  of  the  configuration  vector  fi ,  which  plays  the  role  of  file  force  field. 
Under  typical  conditions  of  TSD  measurement,  when  a  thin  ice  specimen  is  placed  between 
shorted  electrodes  (the  resistance  of  the  measuring  circuit  is  much  less  than  the  resistance  of  ice), 
the  electric  field  within  the  ice  E = 0,  and  the  major  force  driving  chargecarriers  is  Q.  Therefore, 
it  is  more  correct  to  speak  about  relaxation  of  the  configuration  vector  Q  rather  than  relaxation 


of  polarization  P.  Since  we  have 

^  =  ih  -h)  -  ih  -ji)  =  ^ 

dt  ei  63 

(193) 

/ i  =ei{ji  -ji) ;  /b  =  esiis -74) 

(194) 

/i  =  ;  E  =  0;  gradnj  =  0 

(195) 
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a.  Ice  doped  with  HCL  in  different  concentrations. 


b.  Pure  ice  and  ice  doped  with  HCl  and  NH4  OH. 


Figure43.  TSDcurrentsoficespecimens  {after  ZutretskHetal. 

1987b). 

then  instead  of  eq  190  we  obtain  for  current  density 

/=/k,n+7B=4>Q(?-M-  (196) 

\e^  e\  / 

But  diis  is  an  absolutely  amazing  result!  Equatirm  196  predicts  that  in  the  relaxation  process 
dvuing  heating,  the  TSD  current  will  alter  the  direction  at  the  "crossover"  point,  Le.,  where 
T^hat  iS/  relaxation  currents,  initially  flowing  in  the  directicm  opposite  to  the 
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applied  external  field,  after  that  point,  change  direction  and  flow  wifti  it!  This  peculiarity  was 
theoretically  pointed  out  for  the  first  time  by  Zaretskii  et  al.  (1987b),  and  experimentally  the 
direction  change  of  TSD  currents  in  ice  was  reported  by  Johari  and  Jones  (1975).  Later,  this  phe¬ 
nomenon  was  studied  in  detail  by  Aziev  et  al.  (1987).  In  the  report  of  Zaretskii  et  al.  (1987b),  it 
was  also  demonstrated  that  a  similar  current  inversion  can  be  observed  within  another  regime 
of  so-called  Thermally  Stimulated  Currents  (TSC).  Within  this  regime  a  specimen  is  heated  in 
a  strong,  constant  electric  field  and  currents  flowing  dirough  it  are  measured.  In  fact,  such  cur¬ 
rent  inversion  had  already  been  observed  in  the  work  of  Dengel  et  al.  (1964),  in  which  the  alter¬ 
nations  of  current  direction  occurred  at  T  =  112  K  and  T  =  144  K.  At  that  time  such  unusual 
behavior  of  TSC  was  interpreted  in  terms  of  the  transition  of  ice  into  a  ferroelectric  state!  Now, 
however,  in  the  wake  of  numerous  structural  studies,  we  know  that  there  is  no  ferroelectrical 
trarvsition  in  ice  in  ^is  temperature  range.  It  is  easy  to  show,  considering  relaxation  of  rather 
than  relaxation  of  P,  that 

dU  =  -  H  (197) 

dt  XD 

and  instead  of  eq  192  employed  for  determining  of  in  the  TSD  method,  an  essentially  more 

complicated  equation  arises 


1(0-  p(0 

oo 

j/(tj  fiitjdr 


where 


(198) 


^ion  ^  35. 

p  =  Lfi - lii.  (199) 

(26  _  2i2a) 

U3  ei  I 

The  extensive,  accumulated  experimental  data  should  be  processed  using  formula  198  rather 
than  192.  However,  this  rather  time<onsuming  procedure  has  never  been  used.  Still,  if  there  is 
confidence  (from  independent  experiments)  that  the  "crossover"  does  not  occur  within  the  in¬ 
vestigated  temperature  region,  then  we  can  also  use  (with  low  accuracy)  the  more  simple  rela¬ 
tionship  (eq  192). 


CONCLUSION 

At  this  point  we  finish  the  consideration  of  the  electrical  properties  of  ice.  Some  particular 
aspects  that  are  usually  cortsidered  along  with  the  electrical  properties  of  ice  will  be  transferred 
to  futiue  reports  to  avoid  undue  complexity  and  length  of  this  accoimt.  Special  questions  dis¬ 
cussed  herein  are,  as  a  rule,  on  a  frontier  between  various  fields  of  physics.  Thus,  we  will  con¬ 
sider  in  detail  the  surface  conductivity  of  ice  in  the  report.  Surface  of  Ice.  There  the  results  of 
measurements  of  the  Hall-effect  will  also  be  discussed,  since  some  fairly  reliable  data  are  obvi¬ 
ously  related  to  surface  conductivity  only.  Thermoelectric  effects  will  be  coirsidered  in  the  re¬ 
port,  Thentml  Properiies,  along  with  diffusion  processes  in  ice. 

In  the  report  concerned  with  electromechaxucal  effects,  we  will  also  consider  the  effect  of 
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pressure  on  electrical  properties,  pseudo  piezoelectric  properties  and  electromagnetic  emission 
from  cracks  in  ice.  All  these  phenomena  are  interrelated.  In  the  report.  Electro-optica!  Effects  in 
Ice,  both  intrinsic  and  impurity-caused  photoconductivity  of  ice  will  be  described  as  well  as 
photoelectrical  phenomena  at  ice/ electron  conductor  interfaces. 

As  a  unifying  framework  in  the  descriphon  and  discussion  of  the  electrical  properties  of  ice, 
the  author  employed  Jaccard's  theory.  This  seems  to  have  helped  preserve  a  single,  coherent 
framework  in  presenting  the  material.  Such  intensive  use  of  this  theory,  in  spite  of  its  evident 
merits,  does  not  imply  that  Jaccard's  theory  has  no  drawbacks  and  is  able  to  cover  all  "white 
spots"  on  the  ice  physics  map.  Nor  is  it  a  unique  theory.  Jaccard's  theory  is  empirical  or  half- 
empirical,  i.e.,  it  was  not  derived  from  "first  principles."  This  theory  is  based  upon,  firstly,  the 
concept  of  two  fundamentaUy  different  types  of  charge  carriers  (ions  and  Bjerrum-defects)  and, 
secondly,  concepts  about  changes  occurring  in  ice  structure  as  a  result  of  motion  of  these  de¬ 
fects. 

As  far  as  the  first  aspect  is  concerned,  then,  our  ideas  about  the  real  microscopic  structure  of 
ice  defects  may  change  strongly — and  already  are  changing.  The  author  refers  the  reader  to  a 
report  from  the  present  series  called  The  Structure  of  Ordinary  Ice  //, — Part  II.  Real  Structure  of  Ice, 
Defects,  in  which  a  descriphon  is  given  of  alternative  models  of  protonic  structure  defects:  a 
soliton  model,  a  lattice  reconstruchon  model  and  a  model  of  radicals. 

As  far  as  the  second  whole  set  of  concepts  on  which  Jaccard's  model  rests  is  concerned,  I 
believe  it  will  remain  whole,  and  this  make  the  model  extremely  viable.  In  fact,  for  the  validity 
of  the  model,  it  is  important  that  the  hydrogen  bonds  along  which  the  defects  move  would  re¬ 
orient  in  the  manner  perceived  by  Jaccard.  This  will  be  retained  in  all  other  models  (the  soliton 
model,  for  instance).  Thus  the  whole  mathemadcal  descriphon  suggested  by  Jaccard  will  remain! 

The  model  appears  to  be  able  to  account  not  only  for  ac  and  dc  conduchvity  and  dielectric 
permittivity,  but  for  numerous  physical  phenomena  that  appeared  after  the  reports  of  Jaccard 
(1959, 1964). 

This  theory,  of  course,  owing  to  its  empirical  and  schemahc  nature,  cannot  give  answers  to 
all  queshons,  such  as  why  the  mobilihes  of  ions  and  Bjerrum  defects  are  what  they  are.  The 
answer  to  this  and  other  queshons  must  be  obtained  in  the  context  of  theory  shll  in  develop¬ 
ment,  based  upon  the  quantum  theory  of  solids. 
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